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ABSTRACT 

The  purpose  of  this  thesis  was  to  evaluate  the  accuracy  of  the  existing  Mini- 
Ranger  network  in  Monterey  Bay  and  to  suggest  means  by  which  the  network 
can  be  improved.  This  network  consists  of  six  stations  located  around  the  Bay 
and  installed  by  the  Monterey  Bay  Aquarium  Institute  (MBARI)  and  the  Naval 
Postgraduate  School  (NPS). 

In  order  to  undertake  this  task,  data  were  made  available  by  MBARI  which 
were  collected  on  a  cruise  made  by  their  vessel  "Ft.  Lobos".  Additional  data  were 
gathered  on  a  second  cruise  by  the  vessel  "Pt.Sur"  of  the  NPS. 

The  data  analysis  indicated  network  problems.  An  effort  was  made  to  iden- 
tify these  problems  and  to  compute  various  correctors.  In  addition,  an  equation 
has  been  derixed  which  enn'^'es  use  of  Mini- Ranger  data  collected  when  signal 
strengths  are  low. 

The  estimation  of  the  accuracies  obtained  from  the  network  through  the  var- 
ious tests  applied,  and  the  conclusions  drawn,  can  be  used  as  a  guide  to  future 
users  of  the  svstem. 
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I.     INTRODUCTION 

In  1988,  the  Monterey  Bay  Aquarium  Research  Institute  (MBARI)  and  the 
Naval  Postgraduate  School  (NPS)  began  work  to  create  a  permanent  network  of 
Mini-Ranger  stations  around  Monterey  Bay.  The  net  was  to  provide  precise  po- 
sitioning for  research  vessels  and  other  scientific  operations.  Initially  MBARI 
was  seeking  coverage  of  as  much  of  the  bay  as  possible,  especially  the  Monterey 
canyon  feature,  with  absolute  horizontal  position  acuracies  on  the  order  of  a  few 
meters.  It  was  understood  this  would  be  difficult  to  achieve,  given  the  equip- 
ment's specified  range  (40  nm).  This  study  was  undertaken  to  determine  what 
levels  of  accuracy  and  spacial  coverage  were  being  achieved,  as  the  net  was  con- 
figured in  early  1989.  and  also  to  suggest  possible  improvements. 

Positioning  data  from  two  cruises  were  examined,  one  from  a  "Point  Lobos" 
cruise  on  31  MAR  89.  and  a  second  from  "Point  Sur"  cruise  on  22  SEP  89  (the 
"Point  Lobos"  and  "Point  Sur"  are  vessels  used  by  MBARI  and  NPS,  respec- 
tively). Significant  errors  were  found  in  a  high  percentage  of  the  position  fixes 
from  the  "point  Lobos"  data  set.  It  was  apparent  that  at  least  some  of  the  prob- 
lem arose  from  the  use  of  weak  signals  (low  signal  strength  readings)  and  incom- 
pletely calibrated  equipment.  Consequently,  additional  measurements  were  made 
onshore  o\er  known  length  baselines  in  order  to  quantify  the  systematic  and 
random  error  components  of  the  range  measurement  under  \'arious  signal 
strength  conditions. ^Thcsc  findings  stimulate  the  direction  of  the  rest  of  the 
study: 

•  What  were  the  main  sources  of  errors? 

•  Is  there  any  postprocessing  that  can  correct  the  data  and  deri\'e  better  posi- 
tions? 

•  What  is  a  realistic  estimate  of  the  accuracy  achie\'ed? 

•  Are  there  any  chances  for  impro\ement? 


1  It  should  be  noted  that  the  "Point  Lobos"  Master  station  failed  before  the  baseline  testing 
was  run.  Tlius  it  was  not  possible  to  conllrm  the  error  estimates  derived  from  analysis  of  the  cruise 
data  with  independent  baseline  measurements. 


II.     DESCRIPTION 

A.     MINI-RANGER  SYSTEM 

The  Motorola  Mini-Ranger  system  being  used  by  MBARI  and  NPS  is  the 
Falcon  484.  It  is  widely  used  for  precise  positioning  in  a  \'ariety  of  nearshore 
marine  appUcations  like  hydrographic  surveying,  dredging,  oceanographic  data 
collection,  etc.  Numerous  excellent  descriptions  of  the  system  are  a\aliable  in  the 
literature. 

The  system  operates  on  the  principle  of  a  pulse  coded  transponder.  The 
Master  station  (also  refered  as  to  as  the  Receiver-Transmitter)  is  located  on  the 
ship  and  ranges  to  a  maximum  of  four  Reference  stations  ashore.  The  maximum 
update  rate  is  once  per  second.  The  Master  and  Reference  stations  transmit  on 
different  frequencies  in  the  microwaxe  band.  5570  and  5480  Mhz.  respcctixely. 
Both  range  and  signal  strength  data  are  output  for  each  Reference  station  inter- 
rogated. 

Signal  strengths  are  computed  at  the  Master  as  a  relati\"e  indication  of  the 
power  le\el  receixed  from  the  Reference  station.  Values  can  \ary  from  about  5 
to  99.  The  manufacturer  warns  that  any  ranges  with  signal  strength  below  13  are 
"weak"  signals,  subject  to  random  errors  exceeding  the  2  meters  (la)  specifica- 
tion. Signal  strength  ob\  iously  falls  off  with  increasing  range,  but  can  also  drop 
due  to  any  degradation  of  propagation  conditions  between  the  ship  and  shore 
station.  Hence,  weak  signals  can  and  do  occur  at  any  range.  Monitoring  of  the 
recei\"ed  signal  strength  is  critical  to  assessing  the  accuracy  of  all  range  data. 

The  range  is  computed  at  the  Master  station  based  on  the  measured  round- 
trip  travel  time  for  the  pulse  and  a  preset  refracti\e  index  for  microwave  propa- 
gation in  the  atmosphere.  Unless  changed  by  the  user,  the  recei\er  assumes  a 
\alue  of  N  =(n-l)xlO'^  =  320  as  an  indicator  for  the  refractive  index  (n).  which  is 
the  value  used  throughout  this  study.  Output  ranges  also  can  be  corrected  for  any 
constant  (systematic)  error  pre\iously  entered  by  the  user.  Because  this  constant 
error  includes  time  delays  at  both  ends  of  the  line,  it  is  essential  to  ha\c  a  specific 


calibration  of  each  Master/Reference  station  pair.  The  manufacturer  notes  that 
without  cahbration,  systematic  range  errors  of  up  to  ten  meters  are  possible. 

It  is  possible  to  have  the  receiver  correct  ranges  from  "slope  to  horizontal" 
distances  before  output,  based  upon  user  provided  elevations  of  the  Master  and 
Reference  stations.  During  this  study,  both  the  MBARI  and  NPS  equipment 
were  configured  to  output  "slope"  ranges  only,  choosing  to  handle  this  correction 
as  part  of  the  position  computation  algorithm. 

The  Reference  stations  are  supplied  with  electric  current  of  22  to  32  Volts 
DC  from  a  pair  of  common  car  batteries,  or  permanently  from  an  AC  supply 
through  a  converter.  A  possible  insufficient  current  supply  can  cause  errors  in  the 
displayed  distance  in  the  recei\er  [Mini-Ranger  Operation  Manual,  p.II-1 1]. 

The  standard  system  has  a  maximum  range  of  37  km  (20  n.m).  with  Refer- 
ence station  antennas  of  13  db.  The  maximum  range  can  be  extented  to  75  km 
(40  n.m)  by  raising  the  antenna  gain  of  the  Reference  stations  to  19  db,  if  the 
station's  elevation  is  sufficient  (due  to  the  dip  of  the  horizon).  The  master  station, 
which  is  on  the  \essel.  has  an  omni  directional  antenna  of  6  db.  The  discussion 
and  results  in  this  thesis  are  all  gi\en  using  extended  (19  db)  gain  antennas  on  the 
shore  stations. 

Another  significant  factor  for  a  Reference  station  is  the  directi\  ity  of  its  an- 
tenna. The  gi\en  maximum  range  of  the  system  corresponds  to  the  main  lobe  di- 
rection and  drops  as  one  mo\es  to  the  sides.  Figure  1  shows  the  pattern  of  a  19 
db  gain  antenna  of  a  shore  station.  We  can  see  that  the  range  drops  under  30 
km  at  40-  off  the  center  direction.  This  is  significant  for  a  station  that  uses  its 
maximum  range  ability  and  thus  makes  the  choise  o[  antenna  directi\ity  critical. 

B.     MONTEREY  BAY  NETWORK  DESCRIPTION 

The  Motorola  Mini-Ranger  network  in  Monterey  Bay  consists  of  6  stations 
with  high  gain  antennas  (19  db).  installed  permanently  on  known  positions 
around  the  Bay.  The  installation  was  done  by  both  MBARI  and  NPS  ser\ing 
primarily  the  Biological  and  Oceanographic  research  of  their  vessels  "Point 
Lobos"  and  "Point  Sur"  respecti\ely.    The  positions  at  which  the  stations  are  in- 
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Figure  1.        19-db  ANTENNA  PATTERNS. 


Table  1.       PRELIMINARY    POSITIONS    OF    THE    MINI-RANGER 
STATIONS 


Name 

Latitude(N) 

Longitude(W) 

deg  min  sec 

UTM  coordi- 
nates 

Scale 
factor 

Ele- 
vation 
(m) 

Trevor 

36  58  28.785 
122  02  31.536 

585260.161 
4092490.284 

.999690 

51.7 

Wats 

36  54  26.895 
121  50  39.819 

602947.672 
4085231.399 

.999731 

23.1 

Packard 

36  49  44.407 
121  46  04.912 

609863.128 
4076611.345 

.999749 

33.9 

Hays 

36  38  33.826 
121  47  45.895 

607621.289 
4055915.264 

.999743 

137.2 

Doppler 

36  36  05.601 
121  52  37.479 

600434.009 
4051260.014 

.999724 

10.0 

Hank 

36  36  25.381 
121  55  08.722 

596669.490 
4051826.411 

.999715 

134.5 

stalled  are  not  ideal  due  both  to  the  permanent  power  supply  problem  and  to  the 
fact  that  some  of  the  ideal  station  positions  are  on  pri\ate  properties. 

On  board  the  "Point  Lobos"  the  Mini-Ranger  omni-directional  antenna  was 
installed  on  the  mast  9.1  meters  abo\e  the  water  level.  Thus  the  maximum 
achievable  range  due  to  the  line  of  sight  effect  was  increased. 

The  shores  of  the  Bay  and  the  nearshore  areas  are  generally  low  sand  dunes 
so  some  stations  don't  ha\e  the  proper  altitude  to  reach  their  maximum  range. 

Table  1  gives  the  preliminary  positions  [Schnebele,  1989]  of  the  stations,  the 
UTM  grid  coordinates  and  their  scale  factor.  The  horizontal  datum  is  the 
American  Datum  of  1983  (NAD  83).  The  UTM  projection  used  is  for  Zone  10 
(central  meridian  123^  W.)  and  the  scale  factor  0.9996.    Any  errors  in  the  posi- 


tions  of  these  stations  will  propagate  into  the  computed  vessel  position  by  an 
amount  that  depends  on  the  geometry.  If  the  errors  are  smaller  than  se\eral 
decimeters  then  there  is  no  significant  effect.  Revised  positions  are  shown  in  Ta- 
ble 3,  section  IV-C. 

The  stations  are  not  ideally  positioned  and  so  may  be  moved  to  different  po- 
sitions in  the  future  to  allow  better  coverage  of  the  areas  of  interest. 

Station  Doppler  is  not  yet  permanently  installed  for  safety  reasons. 


III.  THEORETICAL  APPROACH  TO  THE  PROBLEM 

A.  THE  RANGING  PROBLEM 

The  technique  used  for  computing  the  ship's  position  using  ranges  from 
known  reference  stations  will  be  derived  and  explained  in  the  beginning  so  the 
reader  will  be  able  to  follow  and  understand  the  data  manipulation  and  results 
of  this  study.  It's  well  known  that  one  way  to  handle  the  computations  invohing 
ranging  data  is  to  work  on  some  appropriate  projection  (e.g  the  UTM  projection). 
There  are  se\eral  techniques  that  can  be  used  for  the  position  derivation: 

•  A\  eraging  the  coordinates  resulting  from  all  possible  combinations  of  meas- 
ured ranges. 

•  Position  calculation  using  the  two  most  accurate  ranges  (using  of  the  rest  as 
a  check). 

•  Graphical  methods. 

In  this  thesis  we  use  the  variation  of  coordinates  method  which  is  described 
in  Cross  [1981].  For  the  position  deri\ation.  at  least  2  ranges  are  needed  (two 
unknown  parameters  in  the  plane  solution  X.Y).  With  more  ranges  one  can  esti- 
mate the  accuracy  of  the  resulting  position  and  also  make  statistical  decisions 
about  the  quality  of  the  data. 

Figure  2  illustrates  the  concept  of  ranging  from  the  i-th  station.  Consider  a 
station  S,  in  a  coordinate  system  X\' .  The  obser\'ed  range,  reduced  to  the  plane. 
is  OR,  (not  shown).  Consider  an  assumed  position  (AP).  The  distance  between 
the  reference  station  and  the  assumed  position  is  CR,.  The  calculation  of  the 
values  dX.dY.  which  represent  the  corrections  that  ha\e  to  be  applied  on  the  AP, 
will  gi\e  the  desired  true  position  (TP)  after  some  iterations. 

As  is  ob\ious  from  the  figure: 

(  ^i)AP  ="  ^  AP~  ^  Si 

The  key  equation  is  the  one  for  the  calculated  range: 


Figure  2.      VARIATION  OF  COORDINATES  METTIOD 


{CR,f  =  {x;)],p  +  {y;)\p 


In  Older  lo  find  Ihe  transformation  from  range  coordinates  to  position  coor- 
dinates it  is  necessary  to  find  the  derivatixe: 

2CR,<iCRi  =  2{X])^ipdX  +2{}))^.^pdY 
dCRi  =  MdX  +  NdY 


where 


The  estimate  of  the  true  range  then  is  equal  to  the  observed  range  plus  an 
error  which  can  be  positi\e  or  negative,  called  the  residual  J',.  It  is  also  equal  to 
the  calculated  range  plus  the  increment  dCR,. 

The  true  range  then  is  OR,  +  I]  where  F,  is  the  residual.  But  it's  also 
CR,  +  dCR,.  Solving  for  dCR„  we  take: 


and  finally: 


dCRi=  {OR  -  CR}^+1) 


M,dX  +  yV,(^}'  =  {OR  -  CR),  +  l) 


For  the  n  ranges  problem,  in  matrix  form,  the  equation  may  be  stated  as  fol- 


lows: 


AX= B+ V 


Where  matrices  A  .B  and  X  respcciixely  are: 


N 


yV, 


5  = 


{OR- 

-CR),' 

{OR- 

-  CR)2 

{OR  - 

-  CR)„ 

x  = 


dX 
dV 


The  standard  solution  for  this  formulation  using  the  least  squares  method 
(Bomford,  1980)  gives: 

X^  =  {a'^IVA}~^A^]VB 


Va  =  AX,  -B 


Where  the  subscript  a  indicates  that  the  value  is  approximate.  The  best  sol- 
ution is  obtained  by  iterating  until  X  converges  to  less  than  some  limit. 

W  is  the  weight  matrix.  The  standard  daviation  o^  of  a  range  measurement 
is  an  expresion  of  its  accuracy.  It  is  standard  to  weight  the  range  measurements 
using  -y  which  ensures  achieving  the  minimum  variance  solution  [Hamilton, 
1964]. 

Assuming  that  there  is  no  correlation  between  measurements,  the  weight 
matrix  becomes  diagonal  as  follows: 


]V  = 


0 
0 


0 

i/cj; 


0 


1^1 


Motorola  estimates  the  standard  deviation  (1  a  )  of  a  single  range  measure- 
ment to  be  2  meters.  The  individual  user  has  to  set  his  own  estimated  (  c  }  for 
his  equipment  (see  section  IV-Hj. 

The  unbiased  estimate  of  the  variance-co\'ariance  matrix  of  the  adjusted  po- 
sition is  gi\en  by 


i:,  =  al{A'^lVAy 


This  matrix  includes  the  variances  and  the  covariances  in  the  coordinate  sys- 
tem we  use  (  X,Y  coordinates)  and  appears  as  follows: 
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^2  <^xy 

\=     J  2  (18) 


Assuming  that  a,  =  Im,  then  the  combination  of  two  ranges  can  result  in  po- 
sitional error  of  7.8  m  (drms)  in  the  limiting  case  of  30°  intersection  angle.  Drms 
is  deri\ed  from  the  unbiased  estimate  of  the  variance  covariance  matrix  as  fol- 
lows: 


\drms  =  -^z  {a^  +  Cy) 

This  error  of  7.8  m  can  be  larger  when  the  system  is  not  calibrated.  For  un- 
cahbrated  systems  the  error  of  a  single  measurement  can  exceed  10  m  [Motorola 
Mini-Ranger  manual]. 

The  quantity  c]  is  the  a  posteriori  variance  of  unit  weight.  It  can  be  computed 
as  follows: 

O^  = 


dof 

Where  dof  is  the  degree  q{  freedom,  or  n-2  when  solving  for  X.Y  position  co- 
ordinates using  n  ranges.  Ideally,  g]  should  equal  unity. 

According  to  Cross  [19S1].  when  the  degrees  of  freedom  are  small  it  is  dan- 
gerous to  use  a  single  determination  <:>{  a^  to  judge  the  quality  o'i  the  position.  It 
is  better  to  use  an  a\'erage  value  of  g^  from  the  whole  data  set. 

But  the  big  question  is  "how  we  can  say  that  the  value  o'i  g^  is  out  of  bounds? 
For  this  question  a  chi  square  one  or  two  tailed  test  is  the  best  indication.  e\en 
though  there  are  situations  where  it  is  not  sufficient  [Uotila.  1975].  If  the  test 
fails,  this  may  be  an  indication  that  one  or  more  of  the  following  error  sources 
exist  that  cause  the  large  value  of  g„. 

•  Error  in  the  mathematical  model  involved  in  the  solution  (Scale  factors, 
Refractixe  Index  used  etc.) 

•  Computational  errors. 
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•  High  correlations  between  observations  or  very  poor  geometry  of  the  situ- 
ation. 

•  Influence  of  the  omitted  higher  order  terms  in  Taylor  series  expansion  which 
can  be  critical  for  higher  order  accuracy  systems  (primarily  when  the  Master 
station  is  far  away  from  to  the  Reference  station). 

•  Incorrect  variance-co\ariance  matrix  of  the  observed  quantities.  It  means 
that  for  a  non  well  calibrated  system,  large  values  of  (j„  are  expected. 

•  Blunders  in  observations  arising  from  various  factors  in  the  measurement 
process  (Multipath,  unstable  readings,  etc). 

So  if  the  statistical  test  used  for  g„  shows  that  something  is  wrong  an  investi- 
gation is  needed  for  the  detection  of  the  possible  error  source  and  its  elimination 
from  the  data. 

B.     PROPAGATION  EFFECTS 

The  line  of  sight  limitation  is  one  of  the  critical  factors  in  the  network  instal- 
lation. The  radius  of  curvature  of^  the  signal  path  is  about  4  times  greater  than  the 
Earth's  curvature,  so  the  heights  of  the  Reference  stations  and  the  height  of  the 
Master  station  antenna  go\ern  the  maximum  range  of  the  system,  according  to 
the  formula.  [Motorola  Mini-Ranger  Manual.  19S1]: 

,7  =  4.04(^7^7  +  ^7^) 

Where  h,  and  /?,„  are  the  heights  in  the  Reference  and  Master  stations  respectively 
and  d  the  maximum  range  in  kilometers. 

The  maximum  range  of  the  system  can  be  reduced  due  to  the  atmospheric 
conditions.  Attenuation  of  the  signal  due  to  rain  is  gi\en  by  [Casey.  1 982]: 

A  —  ar 

Where  A  is  the  attenuation  of  the  signal  in  db,  r  is  the  rainfall  rate  in  mm  hr  and 
a.b  are  frequency  functions.  For  the  Miniranger  mean  frequency  (5.5  GHz),  their 
numerical  \alues  are  a=  1.48  x  10"^  and  b=  1.1469. 

Another  formula  is  gi\'en  for  fog,  based  on  the  water  concentration  per  cubic 
meter: 
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A  =  4.S7xlO~\M/db/km. 


Where  M  is  the  water  content  in  grJnP  and  f  is  the  frequency  in  Ghz. 

Table  2  gives  the  attenuation  of  the  signal  for  rain  and  fog.  at  5.5  Ghz  which 
is  an  average  operating  frequency  for  the  Mini-ranger.  In  terms  of  distance,  we 
have  8.0  km  range  reduction  at  40  km  for  each  db  drop.  So,  for  example,  in  a 
distance  of  40  km  with  heavy  rain  conditions  we  have  1.2  db  drop  in  the  signal 
(1 1.6  km  reduction  in  maximum  range).  This  drop  in  signal  power  can  introduce 
unexpected  range  errors  when  signal  strength  drops  below  the  Critical  Strength 
Threshold  (CST). 

Table  2.       ATTENUATION  OF  THE  SIGNAL  AT  5.5  MHZ 


Precipitation 
Rate  (mm  hr) 

Attenuation 
(db;km). 

Fog  (gr,7;z^) 

Attenuation 
(db/km) 

~) 

0.0015 

1.0 

0.0147 

4  (drizzle) 

0.0072 

1.2 

0.0177 

6 

0.0116 

1.4 

0.0206 

8 

0.0161 

1.6 

0.0236 

10 

0.020S 

1.8 

0.0265 

12 

0.0256 

2.0 

0.0295 

14 

0.0305 

T    -) 

0.0324 

16  (hcaxy) 

0.0356 

2.4 

0.0354 

18 

0.0407 

2.6 

0.0383 

20 

0.0460 

2.8 

0.0413 

A  possible  problem  that  can   appear  in  the  survey  area  is  the  multipath  phe- 
nomenon of  the  signal.    Il  can  cause  constructive  or  destructive  inteferencc  at  the 
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receiver.  There  exist  two  types  of  multipath.  The  first,  vertical  multipath  is 
caused  by  reflection  of  the  signal  on  the  sea  surface.  The  destructive  interference 
at  the  receiver  due  to  vertical  multipath  is  called  "range  hole".  The  second  one 
is  the  horizontal  multipath  due  to  reflection  of  the  signal  on  walls  or  other  sur- 
faces. The  destructive  interference  of  the  horizontal  multipath  is  seen  as  "unstable 
readings". 

Range  holes  are  very  important  in  networks  that  operate  at  line  of  sight  and 
can  be  a  major  source  of  positional  problems.  The  best  way  to  avoid  such  a 
problem  is  to  put  an  additional  Mini-Ranger  Master  station  at  a  different  height 
on  the  \essel. 

The  accuracy  of  the  range  measurements  is  affected  by  the  weather  conditions 
(temperature,  humidity,  pressure).  A  fixed  \alue  of  the  refracti\'e  index  (N)  is  in- 
put into  the  system  prior  to  a  sur\ey  and  kept  constant  throughout  the  survey. 
The  default  refractixe  index  (N  =  320)  corresponds  to  atmospheric  conditions  of 
20-^C  50%  and  1013  mb  respecti\ely. 

Extreme  atmospheric  conditions,  like  high  temperature  along  with  either  high 
relatixe  humidity  or  dry  atmosphere,  can  gi\e  rise  to  large  scale  errors  in  the 
measured  distances.  For  this  reason  further  in\"estigation  was  done  to  estimate 
the  magnitute  of  the  atmospherically  induced  errors  (see  section  IV-E). 

C.     LOW  SIGNAL  STRENGTH 

Se\eral  in\estigators  [e.g  Casey,  1982]  ha\e  shown  that  range  errors,  both 
bias  and  standard  de\iation.  tend  to  increase  once  the  received  signal  strength 
drops  below  some  critical  le\el.  This  le\el,  usually  refered  to  as  the  Critical  Signal 
Threshold  (CST).  is  defined  by  the  users  accuracy  requirements.  Typical  choices 
are  between  signal  strengths  of  8  to  18.  The  Motorola  Falcon  equipment  itself 
flags  signal  strengths  less  than  or  equal  to  13  as  potentially  suspect. 

The  signal  strength  from  the  shore  station  is  also  a  factor  that  can  introduce 
large  errors  in  the  range  measurements,  but  it  is  always  output  by  the  system  and 
can  be  used  in  subsequent  data  processing.  When  signal  strength  is  abo\e  the 
CST  limit,  the  errors  ha\e  a  Gaussian  distribution,  (mean  \alue  zero  and  rela- 
ti\eh'  small  \'ariancc).  But  when  the  signal  strength  is  below  the  CST  limit  the 
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Figure  3.      VARIATION  OF  MEAN  ERROR  WITH  SIGNAL  SI  RENGTII 

errors  follow  a  logarithmic  curve.  That  means  the  range  error  increases  logarith- 
mically as  one  receives  smaller  values  of  signal  strength  [Case\'.  1982]. 

Figure  3  deri\ed  by  Casey,  gives  the  CST  at  SS  =  S  .  For  paiticular  applica- 
tions each  user  can  set  his  own  estimated  \alue  for  CST  in  order  to  either  reject 
the  low  signal  strength  meausurements,  or  to  increase  the  standard  dc\  iation  of 
the  measurement. 

The  SS  of  each  signal  depends  strongly  on  the  distance  from  the  reference 
station,  (it  drops  with  the  in\erse  square  of  the  distance).  It  can  also  \  ary  due  not 
only  to  multipath  phenomenon,  (as  we  explained  earlier),  but  also  due  to  the  an- 
gle to  the  main  lobe  direction  and  the  atmospheric  conditions. 
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IV.     TESTS  ON  THE  MINI-RANGER  NETWORK  IN  MONTEREY  BAY 

A.     MINI-RANGER  CRUISE  DATA 

On  31  March  1989  the  Vessel  Point  Lobos,  from  MBARI,  sailed  for  a  trip  in 
Monterey  Bay  to  make  trials  and  calibration  on  the  Mini-Ranger  stations.  The 
Master  station  was  installed  on  the  mast  at  height  10  m  above  the  water.  During 
this  trip  useful  data  were  collected  and  position  estimates  calculated.  This  data 
was  processed  in  a  preliminary  fashion  in  order  to  detect  possible  error  sources. 

Figure  4  shows  the  routes  of  the  vessel  followed.  A  di\  ision  of  the  Bay  into 
smaller  areas  was  used,  (as  shown  on  the  figure),  in  order  to  subdixide  the  data 
for  processing  and  error  analysis  purposes.  The  division  was  based  on  the  con- 
cept of  nearshore  and  open  b:  y  regions. 

Another  trip  was  subsequently  done  for  NPS  research  purposes  on  22  Sep- 
tember 1989,  and  some  more  data  collected  as  shown  on  the  same  figure.  The 
master  station  for  this  cruise  was  different  and  the  data  w  collected  under  differ- 
ent conditions  of  atmosphere,  time  and  instrumentation.  It  was  considered  possi- 
ble that  the  positional  results  from  these  two  data  sets  may  \ary  in  accuracy. 

First  we  examined  the  earlier  data  set  of  MBARI,  which  already  had  the 
positions  and  their  statistics  calculated.  We  were  informed  by  MBARI  that  for 
the  position  deri\ation  they  used  the  method  o(  \ariation  of  coordinates.  No  de- 
tails on  the  specific  code  in  this  software  were  available. 

The  first  step  was  a  general  \'ic\v  on  the  residuals  and  the  "standard  de\1- 
ation"  provided  in  MBARIs  data  of  the  whole  trip.  This  showed  that  something 
did  not  work  well  when  the  vessel  was  nearshore  ,  especially  in  the  SW  portion 
of  the  Bay.  Figure  5  shows  the  areas  of  the  greater  and  more  \ariable  standard 
dex  iation  of  the  data  set.  In  general  as  shown  in  Figure  5,  areas  E  and  G  ha\"e 
fewer  problems.  It  was  not  ob\ious  which  station  or  stations  combination  was 
responsible  for  the  big  residuals. 

Another  problem  was  the  variation  of  signal  strength  and  its  unexpectedly 
high  \'alues.  In  many  portions  of  the  data  set,  the  signal  strength  was  much  higher 
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Figure  4.      SHIP  ROUTES  FOR  THE  TWO  DATA  SETS/AREA  DIMSION 
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Figure  5.      AREAS  OF  GREATER  STANDARD  DEVIATION 
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than  expected  for  such  a  range.  The  appearance  of  signal  strengths  greater  than 
70  in  distances  more  than  10000  meters  is  unusual.  This  is  a  strong  indication  of 
possible  S3"stem  malfunction. 

Because  the  data  set  did  not  show  the  expected  accuracy,  the  following  series 
of  tests  and  checks  were  performed: 

•  A  visit  to  the  stations  for  inspection  of  directional  limitations,  antenna  in- 
stallation problems  and  other  possible  affecting  factors. 

•  A  check  of  the  station  positions  with  geodetic  survey  techniques,  (triangu- 
lation.  traverse,  or  differential  GPS),  and  a  least  squares  adjustment  of  the 
sur\"ey  data  with  fixed  geodetic  control  points  around  the  Bay. 

•  Check  of  the  software  used  for  the  position  derivation  using  the  variation  of 
coordinates  algorithm  with  accurate  grid  scale  and  slope  distance  correction 
techniques. 

•  A  theoretical  view  of  the  geometrical  accuracies  with  the  use  of  various  sta- 
tion combinations  in  the  Bay  area. 

•  A  check  on  the  effects  of  extreme  meteorological  conditions  to  check  for  the 
possibility  of  introducing  range  errors  in  excess  of  1  meter. 


•  A  range  hole  prediction  for  the  network. 

•  An  in\"esiigation  into  the  possibility  of  correcting  those  ranges  with  signal 
strengths  below  the  CST. 

•  A  calibration  test  for  determining  biased  or  linearly  varying  errors  and  their 
\ariability  with  signal  strength  above  the  CST. 

B.     INSTALLATION  CHECK  AND  ANTENNA  DIRECTIMTIES 

On  20th  of  July,  a  visit  to  the  stations  TREVOR.  WATS  and  PACKARD 
was  done  for  \isual  examination,  check  on  the  antenna  directi\ity  .  possible  sig- 
nal blocking  from  nearby  obstructions  and  other  possible  affecting  factors.  The 
stations  HA'l'S  and  HANK  were  visited  on  the  1st  of  August  and  25th  of  Sep- 
tember respecti\"ely.  The  alignments  of  the  antennas  were  determined  to  approx- 
imately 2^'  using  a  small  portable  magnetic  compass.  All  the  magnetic  bearings 
are  gWen  as  true  bearings.  A  short  description  is  gi\'en  below  for  each  station  and 
a  map  is  presented  (Figure  6)  for  clarification: 

•  TRE\'OR:  This  station  is  installed  on  the  roof  of  a  house  in  Santa  Cruz.  The 
antenna  is  attached  to  the  chimney  of  the  house  and  is  oriented  towards 
175'.  The  signal  is  blocked  by  a  rocky  hill  from  215'-'  to  approximately  100^ 
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Figure  6.        OBSERVED    ANTENNA    DIRECTIVITIES    OF    MINI-RANGER 
NETWORK  IN  MONTEREY  BAY 
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(clockwise).  There  is  no  clear  view  of  the  Northern  near  shore  Bay  region  but 
this  area  is  of  little  interest. 

•  PACKARD:  Installed  on  the  roof  of  a  wood  craft  building,  hooked  on  a  thin 
mast,  on  private  property  between  Watsonxille  and  Moss  Landing.  Its  an- 
tenna points  to  250"  and  may  be  blocked  by  the  Moss  Landing  electric  plant 
stacks  between  202"  and  205"  There  were  no  other  problems  observed  for  this 
station. 

•  WATS:  Installed  on  the  roof  of  a  commercial  building  near  the  beach  at 
Watson\ille.  Its  antenna  points  at  215"  .The  entire  Bay  area  is  unobstructed. 
The  antenna  is  approximately  2  feet  from  the  vertical  wall  of  the  building. 

•  HA'l'S:  Installed  on  the  roof  of  Silas  B.  Hays  hospital  at  Fort  Ord  with  an 
uninterupted  Bay  view  and  directivity  of  its  antenna  to  295"  .  The  perfect 
view  along  with  the  height  of  the  station  makes  it  useful  in  some  portions 
SW  of  Point  Pinos  .  The  signal  is  blocked  from  the  higher  Pacific  Gro\e  hills 
from  253"  back  to  the  South. 

•  HANK:  Installed  on  the  roof  of  a  house  in  Pacific  Gro\e.  While  it  has  good 
height  the  signal  cannot  reach  the  nearshorc  portion  from  Point  Pinos  to 
Monterey  harbour  due  to  high  trees  existing  downhill.  The  directi\"ity  of  its 
antenna  is  345"  and  the  signal  may  been  blocked  by  high  trees  from  005"  to 
035".  The  station  is  offset  from  the  position  specified  by  the  preliminary  co- 
ordinates. This  particular  correction  is  given  in  the  discussion  in  section  C 
below. 

•  DOPPLER:  This  station  was  installed  near  NPS  Beach  Lab.  but  due  to 
safety  reasons  was  subsequenly  removed. 

It  was  not  possible  to  check  the  electric  current  supply  of  the  stations,  but  a 
voltage  stabilizer  can  help  the  system  to  be  kept  healthier  and  unstable  outputs 
can  be  a\oided  in  the  future. 

C.     POSITION  CHECK 

Geographic  coordinates  for  the  station  sites  had  been  determined  on  the 
North  American  Datum  of  1983  (NAD  83)  by  geodetic  sur\ey  techniques,  with 
ties  to  existing  horizontal  control  points  published  by  the  National  Geodetic  Sur- 
vey [Schnebele,  1989].  Conventional  ground  sur\ey  and  Global  Positioning  Sys- 
tem (GPS)  techniques  were  both  used.  The  sur\'ey  data  had  been  analyzed  with 
a  three-dimentional  least-squares  adjustment  software  package  (GEOLAB,  by 
GEOSurw  Inc.).  Absolute  horizontal  position  accuracies  estimated  by  the  ad- 
justment were  better  than  8  cm  (circular  error.  95^''()  confidence  k\e\).  Additional 
work  to  confirm  the  offset  positions  of  the  Mini-Ranger  units  at  PACKARD  and 


WATS,  has  yet  to  be  completed,  but  is  not  expected  to  shift  their  preliminary 
positions  by  more  than  a  few  centimeters  at  most. 

These  station  coordinates  were  converted  to  Universal  Transverse  Mercator 
grid  coordinates  (UTM  Zone  10)  for  use  in  the  vessel  position  computation  algo- 
rithms. The  UTM  system  has  been  adopted  as  the  standard  grid  used  by  MBARI 
and  NPS  with  this  network.  The  specific  Transverse  Mercator  Projection  algo- 
rithm [Floyd,  1985]  is  accurate  to  within  1  cm  for  the  entire  zone. 

As  described  earlier,  station  HANK  was  offset  from  the  preliminary  position. 
The  offset  was  measured  by  a  simple  magnetic  compass  for  azimuth  and  a  tape 
for  distance  and  should  have  errors  no  greater  than  5  cm.  The  measured  bearing 
was  066''  true  and  the  distance  41  inches.  The  two  positions  were  at  the  same 
height,  so  the  corrections  in  UTM  are:  dx  =  0.952,  dy  =  0.424,  dz  =  0.00  in  meters. 

The  results  are  sumarized  and  the  revised  station  positions  provided  in  the 
Table  3.  With  the  exception  of  station  HANK,  the  revised  positions  of  the 
stations  differ  by  less  than  10  cm  from  the  preliminary  positions.  It  is  concluded 
that  the  station  positions  are  not  the  source  of  the  observed  inaccuracies. 

D.     ESTIMATION  AND  COMPARISON  OF  DERIVED  POSITIONS 
1.     Development  of  position  derivation  algoiithni 

The  measured  ranges  were  used  to  recompute  the  ship's  positions  by  the 
method  of  \ariation  of  coordinates  [Cross,  1981].  These  recomputed  positions 
were  then  compared  with  those  pro\ided  by  MBARI.  The  comparisons  showed 
that  there  is  a  discrepancy  between  the  given  positions  and  the  recomputed  posi- 
tions on  the  first  data  set.  which  in  some  cases  were  significant.  Because  of  these 
discrepancies,  we  ignored  all  the  gi\en  MBARI  positions  and  recomputed  all  their 
data. 

For  this  purpose,  a  program  was  written  in  Turbo  Basic.  This  program 
uses  the  \ariation  of  coordinates  method  on  the  UTM  projection  surface.  Scale 
factors  were  computed  for  each  line.  The  program  is  gi\en  in  the  Appendix  I  and 
can  be  followed  by  comments  pro\ided  in  the  code. 

The  measured  distances  from  the  system,  for  each  case,  are  corrected  for 
their  slope  with  the  assumption  of  a  right  triangle  (ignoring  the  cur\aiure  of  the 
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Table  3.       REVISED         POSITIONS 
STATIONS 


OF 


MINI-RANGER 


Name 

Latitude(N) 

Longitude(W) 

deg  min  sec 

UTM  coordinates 

Elevation 

Trevor 

36  58  28.785 
122  02  31.536 

585260.158 
4092490.287 

51.1 

Wats 

36  54  26.896 
121  50  39.821 

602947.636 
4085231.429 

22.6 

Packard 

36  49  44.404 
121  46  04,911 

609863.147 
4076611.241 

33.4 

Hays 

36  38  33.82- 
121  47  45.894 

607621.298 
4055915.205 

136.7 

Hank 

36  36  25.393 
121  55  08.683 

596670.450 
4051826.796 

134.1 

Earth).  This  assumption  inserts  an  error  of  4  to  5  cm  at  the  extreme  ranges  which 
is  insignificant  for  Mini-ranger  measurements  in  the  areas  co\ered  b}"  the  net- 
work, but  it  makes  the  calculations  simpler.  For  the  complete  correction  one  can 
look  in  chapter  15  of  Laurila  [1981]. 

The  next  correction  that  was  made  to  the  measured  distances  was 
for  the  scale  factor  of  the  UTM  projection.  The  following  formula  computes  the 
scale  factor  for  the  line  between  the  Master  and  each  Reference  station: 


m  =  A77^(l   + 


6R- 
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Where  m  is  the  required  scale  factor,  m^  is  the  scale  at  the  central  meridian 
(0.9996),  R  is  the  mean  radius  of  the  Earth  (6375  km)  and  £^  is  the  mean  Easting 
value  for  each  case.  The  mean  Easting  value  is  derived  from  the  formula: 

Where  £,  and  E^  are  the  true  Eastings  {x,  -  500000),  of  the  Master  and  the  Ref- 
erence stations  respectively.  The  formula  is  a  simplification  of  the  one  given  by 
Bomford  [1980,  pi 93]  and  the  error  due  to  this  is  of  millimeter  order,  which  for 
our  purposes  is  negligible. 

The  outputs  of  the  program  are  the  least  squares  estimate  for  the  position 
of  the  vessel,  the  variance  covariance  matrix  L^,  the  residual  \ector  (V)  and  the 
a  posteriori  variance  of  unit  weight  ol.  Measurements  were  assumed  independent 
so  the  diagonal  weight  matrix  is  treated  as  a  vector  in  the  calculations. 

The  a  posteriori  \ariance  of  unit  weight  0%  is  the  one  that  is  used  for  as- 
sessing both  the  accuracy  of  the  network  and  the  ranging  data.  0%  can  be  scaled 
as  is  well  known,  by  increasing  the  standard  deviation  of  a  single  range  meas- 
urement. The  resulting  c\  is  examined  statistically  with  the  chi  squared,  two 
tailed  test,  at  a  95-'o  confidence  le\el.  The  null  hypothesis.  gI=  1,  was  tested 
against  the  alternative  cli^  1.  When  the  null  hypothesis  is  accepted,  it  is  an  in- 
dication, assuming  no  data  blunders,  that  the  system  works  within  the  specifica- 
tions pro\'ided  by  the  manufacturers.  When  the  null  hypothesis  is  rejected  se\eral 
possible  reasons  for  this  exist.  These  were  summarized  in  section  Ill-A. 

It  also  should  be  noted,  however,  that  the  presence  of  any  unresohed 
systematic  effects  tends  to  in\  alidate  this  type  of  test.  One  potential  source,  not 
modelled  in  the  position  computation  algorithm,  in\"ohes  the  non-simultaneity  of 
range  measurements  on  a  mo\'ing  ship.  In  effect,  the  antenna  is  not  stationary 
over  a  cycle  of  range  measurements  to  all  Reference  stations.  Gi\"en  each  range 
measurement  takes  about  55  milliseconds,  four  range  measurements  may  be 
spread  o\"er  a  period  of  about  0.15  seconds.  At  a  nominal  \essel  speed  of  6  knots, 
the  possible  error  is  about  0.45  cm.  which  is  of  little  importance  in  this  analysis. 
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Additional  error  can  be  inserted  due  to  mast  motion.  In  rough  sea  condi- 
tions the  error  is  big  enough  (more  than  1  meter)  and  it  has  to  be  taken  into  ac- 
count by  increasing  the  estimated  standard  deviation  of  the  range  measurements. 
2.     First  data  set 

The  results  of  the  chi  square  test  at  a  95%  confidence  level  using  standard 
deviations  of  2  and  3  meters  on  the  measured  ranges  respectively,  appear  in  Table 
4. 

The  acceptable  values  of  ol  under  these  conditions  are: 

For  three  ranges  0.00  <  ul  <  5.024 

For  four  ranges   0.025  <  u;  <  3.69 

The  numbers  gi\en  in  the  first  two  columns  of  Table  4  represent  the  per- 
centage of  the  number  of  positions  for  which  the  null  hypothesis  was  rejected 
against  the  total  number  of  the  positions  examined.  The  number  pro\idcd  does 
not  include  some  abnormaly  high  a  posteriori  variances  of  unit  weight  which  were 
caused  by  multipath  effects,  or  low  SS.  These  data  were  rejected  as  blunders. 
Their  percentage  appears  in  the  last  column. 

As  can  be  seen.  e\en  when  a  =3m.  the  number  of  rejections  is  clearly 
much  higher  than  the  5"o   which  one  would  normally  expect. 

In  the  column  'Difference'  we  can  see  how  many  of  the  positions  that 
were  rejected  when  a  =2m.  are  accepted  with  a  =3m.  This  is  an  indication  ei- 
ther of  poor  calibration  of  the  system  or  that  the  system  does  not  perform  to  the 
standard  claimed  by  the  manufacturer.  It  should  be  noted  that  the  range  data 
had  been  corrected  for  an  estimate  of  systematic  errors  determined  in  a  prior 
calibration. 

The  results  shown  in  Table  4  suggest  the  possibility  that  the  master  sta- 
tion used  may  ha\e  beeen  poorly  calibrated.  For  this  reason,  an  examination  of 
the  distribution  of  the  residuals  for  each  reference  station  was  undertaken.  After 
examining  the  distribution  of  the  range  residuals  for  each  reference  station,  it 
became  clear  that  some  stations  were  consistently  biasing  the  posiuon  solutions. 
From  the  distribution  of  these  residuals,  the  bias  was  estimated  by  comparing 
their  mean,  median  and  most  probable  \alues.     The  mean  residuals  and  the 
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Table  4.       REJECTED     PERCENTAGE     WITHOUT     REVISED 
CORRECTORS 


AREA 

(7=2 

(7    =    3 

Difference 

Rejected 
blunders 

A 

68"  0 

48"  0 

20"/o 

\\% 

B 

82"^o 

37"/o 

45»/o 

9"/o 

E 

69"  0 

22"^o 

4  7  "4 

25"  0 

F 

27"  0 

13"o 

14"o 

33"o 

G 

43"o 

27% 

16"o 

29% 

adopted  correctors  are  shown  in  Table  5.  This  simple  averaging  of  obser\  ed  res- 
iduals is  not  an  ideal  procedure  for  estimating  bias  errors.  The  least-squares 
process  tends  to  redistribute  an  error  in  any  one  range  to  errors  (residuals)  in  all 
ranges.  Recognizing  this  limitation,  howe\er,  the  procedure  gives  a  useful  esti- 
mate of  uncorrected  bias  errors. 

The  results  for  the  corrected  data  are  shown  in  Table  6  and  can  be  di- 
rectly compared  with  Table  4.  The  improvement  as  we  see  is  abo\e  65  ", o. 

In  the  column  'Difference'  we  can  see.  in  comparison  with  the  one  of  table 
4.  that  the  \  alues  are  much  lower,  which  is  an  indication  that  most  of  the  cali- 
bration uncertainty  has  been  remo\ed. 

E\en  after  this  procedure,  some  a  posteriori  \ariances  of  unit  weight  are 
still  high,  especially  in  areas  A  and  B.  The  possible  reasons  for  this  can  be: 


Incomplete  remo\al  of  the  blunders. 

Operation  of  some  stations  at  low  SS  which  in  turn  causes  larger  deviations 
from  the  normal  values. 

Incomplete  remo\al  of  the  bias  errors  described  abo\e.  The  technique  used 
does  not  eliminate  the  entire  error,  but  only  estimates  it. 

For  any  one  of  a  number  of  reasons  the  signal  can  be  lost  temporarily. 
When  it  is  regained  it  sometimes  lacks  stabilitv  for  se\eral  readimzs. 
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Table  5.       MEAN  RESIDUAL  AND  VALUES  USED  FOR 
RANGE  CORRECTORS 


STATIONS 

Mean  residual 

Used  corrector 

Ti"e\or 

3.91  m 

-3.5  m 

Wats 

-1.66  m 

1.5  m 

Packard 

-0.24  m 

0.0  m 

Hays 

2.61  m 

-2.5  m 

Doppler 

3.64  m 

-3.5  m 

Hank 

-0.24  m 

0.0  m 

Table  6.       REJECTED   PERCENTAGE   AND   IMPROVEMENT  WITH   RE- 
VISED CORRECTORS 


AREA 

o  =  2 

u  =  3 

Difference 

Improve- 
ment 

A 

38"  0 

22"o 

16"o 

44"  0 

B 

31  «o 

14"o 

17% 

62  "  0 

E 

14"o 

06"  o 

8"o 

79  "o 

F 

8% 

0"'o 

8% 

70"  0 

G 

13"o 

7"o 

6"o 

69"  0 

•    Higlier  \alue  of  tlie  standard  deviation  of  tlic  system  used  than  the  one  used 
for  the  weight  matrix  formation. 
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3.     Second  data  set 

The  second  data  set  was  examined  statistically  in  the  same  way  and  the 
results  appear  in  Table  7.  We  can  see  that  the  second  data  set  gives  much  better 
solutions.  This  was  expected  because  it  was  known  that  the  master  station  used 
for  this  data  collection,  was  well  calibrated. 

Table  7.       REJECTED  PERCENTAGE  FOR  THE  SECOND 
DATA  SET 


AREA 

Rejected  percentage  for  cr  =  2 

A 

9% 

B 

(V^o 

C 

6% 

D 

7"o 

E 

7^) 

F 

7% 

The  pre\  iuus  data  set.  for  reasons  noted  earlier,  showed  major  problems 
when  the  ship  was  in  areas  A  and  B.  The  second  data  set,  howe\er,  shows  con- 
sistent results  in  all  areas  and  a  low  rejection  percentage  even  when  c  =2m. 
These  results  suggest  that  the  standard  deviation  of  the  range  measurement 
should  be  between  2  and  3  meters. 

E.     GEOMETRICAL  ACCURACIES 

The  position  quality  given  from  the  Mini-ranger  network  in  Monterey  Bay 
depends  strongly  on  the  relati\  e  geometry  of  the  ship  and  the  reference  stations. 
The  geometric  position  quality  can  be  represented  in  terms  of  error  ellipses,  con- 
fidence circles,  or  with  se\eral  other  appropriate  techniques.  When  se\eral 
stations  are  used,  computer  techniques  are  necessary  if  detailed  analysis  is  to  be 
undertaken.  Such  techniques  are  described  in  Perugini  [1988]. 
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The  provided  figures  in  Appendix  2  were  taken  from  computer  calculations 
using  software  associated  with  MBARI's  Navigation  processor.  The  numbers 
provided  represent  the  semi-major  axis  of  the  predicted  error  ellipse  at  95%  con- 
fidence level  that  is  formed  at  the  particular  position,  assuming  the  standard  de- 
viation of  a  single  range  measurement  is  3  meters. 

The  best  working  combination  appears  in  Figure  7  and  is  obtained  using 
stations  TREVOR.  PACKARD,  HAYES  and  HANK.  We  can  see  that  the  ge- 
ometry of  the  position  is  very  good  with  two  combinations  of  only  three  stations, 
i.e;  station  combinations  TREVOR,  PACKARD,  HANK  and  TREVOR, 
PACKARD,  HAYES  (see  Appendix  2) 

As  is  easily  seen,  station  TREVOR  occupies  a  position  that  gi\es  strong  ge- 
ometry in  the  mid-Bay  regions  and  its  position  is  the  most  critical  for  the  network. 
On  the  North  side  of  Monterey  Bay  there  are  also  other  higher  positions  on  which 
a  permanent  station  could  be  installed  with  a  much  better  view  of  the  Bay.  Most 
of  them  are  in  the  Uni\ersity  of  California  Santa  Cruz  campus  and  some 
arrangments  ha\e  to  be  made. 

In  the  \Vatson\ille  area  the  hills  over  Buena  Vista  Point  ha\e  a  good  Bay 
\iew  but  the  problem  of  a  permanent  electrical  supply  and  security  has  to  be 
sohed.  Such  a  position  would  o\"crcome  the  height  problem  of  the  WATS  station 
while  maintaining  its  geometry. 

Stations  PACKARD  and  HAYS  are  at  the  highest  elevations  in  their  areas. 
There  is  no  adxantage  to  shifting  their  positions  without  deteriorating  the  geom- 
etry. 

Station  DOPPLER  does  not  add  much  to  the  geometric  quality  of  the  net- 
work in  the  mid-Bay  but  can  be  used  for  positioning  near  Monterey  harbor  where 
signals  from  PACKARD  and  HANK  are  obstructed.  It  also  could  be  shifted  NW 
of  Point  Pinos  on  the  lighthouse  or  on  another  high  position  to  impro\e  the  ge- 
ometry in  the  Monterey  canyon  area.  Further  insestigation  on  this  station  has 
to  be  done  for  the  best  choice  of  a  permanent  position. 
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HAYS 


Monterey  Bay  Vicinity 


Figure  7.        BEST    WORKING    COMBINATION    OF    THE    MINI-RANGER 
STATIONS 
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Station  HANK  gives  very  good  geometry  but  its  position  has  some  limitations 
because  of  obstructions  to  its  Bay  view.  Some  station  at  the  southern  site  of  the 
Bay  is  necessary  to  maintain  a  good  network  geometry  for  the  areas  of  interest. 

Generally  we  can  say  that  the  geometry  of  the  network  is  not  responsible  for 
the  big  standard  deviations  observed  in  most  regions  around  the  Bay,  but  it  be- 
comes more  important  as  one  mo\es  offshore  near  the  Northen  and  Southern  re- 
gions. 

F.     ATMOSPHERIC  CONDITION  EFFECT  OF  THE  AREA  ON  THE 
MICROWAVE  PROPAGATION  VELOCITY 

An  in\estigalion  was  carried  out  to  see  the  maximum  possible  error  due  to 
propagation  \elocity  which  could  occur  under  extreme  atmospheric  Bay  condi- 
tions as  compared  with  mean  Bay  conditions.  The  atmospheric  formulation  is 
described  in  Chapter  6  p.  133  oi  Laurila  [1981].  The  effects  of  temperature  (T). 
pressure  (P)  and  humidity  (E)  on  the  refracti\e  index  (N)  are  gi\en  from  the 
Essen  formula: 

A'  =  77.62  Y  -  (  -^-^Y^  -  37.19  ^  )E 

Where  P  and  E  arc  in  mb  and  T  is  the  dry  bulb  temperature  in  K'\  But 

iRH)E. 


E 


100 


Where  RH  is  the  relati\e  humidity  (%)  and  E^  is  the  saturation  \apor  prcssuie, 
gi\en  in  mb.  at  the  dry  bulb  temperature.  E^  is  a  logarithmic  function  of  temper- 
ature so  the  effect  on  the  refracti\e  index  depends  on  the  temperature  and  the 
relati\"e  humidity  as  shown  in  Figure  8,  for  various  relative  humidities  and  con- 
stant atmospheric  pressure  at  1013  mb. 

Daily  weather  logs,  kept  by  the  Meteorology  Department  of  the  NPS.  were 
searched  for  the  maxima  and  minima  o'i  temperature  and  relative  humidity  in 
Monterey  during  1988-89.  This  was  done  so  as  to  compute  the  a\erage  \'ariability 
of  the  refractiN  e  index  and  the  possible  error  in  distance  that  is  introduced  for  the 
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Temperature  (C) 


Fiauie  8.      VARIATION  OF  THE  REFRACTIVE  INDEX 


actual  atmospheric  conditions.  The  base  for  the  comparison  is  the  refractive  index 
\alue  (N  =  320)  that  is  used  from  the  Mini-ranger  for  the  range  calculation. 

In  Monterey  Bay  during  the  spring,  temperature  a\erages  about  53  'F  and 
humidity  \  aries  from  30%  during  the  hotest  days  to  95%  during  the  cold  nights. 
During  working  hours,  a  relati\e  humidity  of  55%  with  temperatures  70  P  are 
usual.  Extreme  \  alues  of  95%  RH  and  55  "F  are  rare,  but  if  they  are  accompa- 
nied with  high  atmospheric  pressure,  1030  mb,  then  they  gi\e  N  =  343.6  which 
in  turn  corresponds  to  a  distance  error  of  23.6  ppm.  In  a  possible  distance  of 
40000  m  the  introduced  error  is  about  1  meter.  On  the  other  hand,  low  atmo- 
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spheric  pressure  of  1000  mb  with  80  ""F  and  20%  humidity,  gives  N  =  288  and 
32  ppm  error,  which  in  distance  terms  means  1.28  m  at  40000  m  measured  dis- 
tance. A\erage  \alues  of  N  for  the  spring  during  working  hours  are  between  310 
and  330,  which  indicates  that  the  assumed  Refractive  Index  of  320  is  fair. 

During  the  summer  the  increased  temperature  is  accompanied  with  a  drop  in 
the  relative  humidity  so  the  result  tends  to  be  the  same  as  in  spring  with  an  av- 
erage N  =  320  to  325.  The  extreme  cases  here  also  rarely  introduce  errors  more 
than  1  meter  at  distances  of  40000  m. 

The  fall  season  is  the  one  with  a  slightly  different  mean  N  and  more  vari- 
ability. Warm  weather  with  high  humidity  is  more  likely  to  occur.  In  the  extreme 
cases,  at  maximum  distance,  the  error  could  reach  2  meters,  but  it  rarely  occurs. 
The  average  N  =  330  is  better  to  use. 

Winter  shifts  the  refracti\'e  Index  back  to  the  lower  values  due  to  the  lower 
temperatures  of  the  atmosphere,  which  makes  the  extreme  cases  more  rare.  A 
mean  \alue  of  N  =  315  is  better  for  use. 

Generally  the  resulting  error  should  be  less  than  a  meter  at  extreme  ranges 
and  se\eral  decimeters  on  aNcrage,  when  extreme  atmospheric  conditions  happen. 
Even  though  the  error  is  small,  the  use  of  the  seasonal  refracti\"e  index  is  advised. 

G.     RANGE  HOLES 

Range  holes  occur  due  to  the  vertical  multipath.  which  is  caused  from  re- 
flection of  the  signal  from  the  sea  surface.  This  is  common  to  all  microwa\e  sys- 
tems that  use  a  line-of-sight  propagation  path.  Figure  9  [Gilb-Weedon.  1976] 
illustrates  the  situation. 

We  count  the  range  holes  from  the  first  to  the  Nth.  Here  N  represents  the 
integer  number  of  wavelengths  difference  between  the  direct  and  the  reflected 
paths.  The  first  range  hole  occurs  at  the  longer  distance  from  the  shore  station. 

Se\eral  factors  affect  the  range  hole  characteristics: 

•  The  distance  between  the  reference  and  the  master  stations.  The  longer  the 
distance,  the  wider  the  range  hole  zone. 

•  The  height  of  the  stations.  The  higher  the  reference  stations  the  fewer  the 
ranse  holes  near  them. 
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Figure  9.      RANGE  HOLE  OCCURANCE  CONCEPT 


•  The  wavelength  of  the  signal.  This  affects  both  the  width  of  the  range  holes, 
and  the  distance  from  the  reference  station  at  which  they  occur. 

•  The  atmospheric  conditions  along  the  signal  path  cause  loss  of  SS.  Thus  the 
reflected  signal  from  the  sea  surface,  due  to  additional  losses  from  the  re- 
flection, is  too  weak  to  interfere  the  direct  one,  when  they  both  arri\e  to  the 
master  station. 

•  The  sea  surface  roughness.  This  affects  the  type  of  reflection  which  occuis. 
In  a  strongly  wave  modulated  sea,  the  appearance  of  the  range  holes  is  rare. 

•  Changes  in  the  sea  le\el  due  to  tides,  winds  etc.,  cause  range  holes  to  shift 
location  horizontally. 

Small  reflection  angles  cause  more  frequent  range  holes  because  the  signal  has 
lesser  losses  in  comparison  with  the  case  when  the  reflection  angle  is  bigger.  The 
additional  losses  of  SS.  in  the  second  case,  weaken  the  reflected  signal  vvhich  can 
not  then  interfere  with  the  direct  one. 

In  reality,  the  destructive  interference  is  rarely  total  due  to  losses  of  the  re- 
flected signal  from  the  sea  surface,  but  it  can  reduce  the  SS  to  below  the  sensi- 
ti\ it\'  of  the  recei\er  or  below  the  GST  in  the  case  of  Mini-ranger.  This  will  result 
in  signal  loss,  or  worse,  in  a  bigger  standard  de\iation  of  the  measurement. 
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In  the  Monterey  Bay  network,  the  range  hole  predictions  are  shown  in  Table 
8  for  each  station.  For  the  maximum  range  of  each  station  we  assume  a  ship  an- 
tenna of  10  m. 

Range  holes  closer  than  the  third  occur  rarely  because  of  the  bigger  angle  of 
incidence  which  attenuates  more  of  the  reflected  signal. 

Table  8.       RANGE  HOLE  PREDICTION  FOR  MONTEREY  BAY 


STATIONS 

Maxi- 
mum 
range 

Range  for  the 
first  hole 

Range  for  the 
second  hole 

Range  for  the 
third  hole 

Trevor 

42000  m 

18972  m 

9486  m 

6324  m 

Wats 

30000  m 

8477  m 

4238  m 

2826  m 

Packard 

37000  m 

12440  m 

6220  m 

4147  m 

Hays 

61000  m 

50348  m 

25174  m 

16783  m 

Doppler 

26000  m 

3670  m 

1835  m 

1223  m 

Hank 

60000  m 

49358  m 

24679  m 

16453  m 

Table  8  is  based  on  the  following  formula  [Gilb-Weedon.  1976].  in  which  a 
plane  earth  is  assumed.  R  represents  the  distance  from  the  reference  station  to 
the  ran.ce  hole  center  in  meters. 


R  =  2l:,^ 


The  distance  of  the  maximum  expected  signal  is: 


h- 


'   {2k  +  1) 


Where  /z,  and  //;  are  the  heights  in  meters  of  the  reference  and  the  master  stations 
respecti\"cly.  1  is  the  wa\clength  and  k  is  an  integer  (k  =  1.2,3....). 
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Figure   10.      RANGE  HOLE  SHIFT  VS  TIDAL  HEIGHT 

In  the  Bay  there  is  a  sea  level  variation  due  to  tidal  effects  of  6  to  7  feet 
maximum.  This  effect  tends  to  shift  the  range  hole  center  towards  the  shore  sta- 
tion at  high  water  and  away  from  the  shore  station  at  low  water.  The  o\crall  ef- 
fect of  this  phenomenon  for  the  stations  appear  in  Figure  10.  for  each  range  hole. 

Typical  wind  and  wave  conditions  in  the  Bay  are  such  that  'perfect'  signal 
reflections  are  not  expected.  Thus  range  hole  problems  should  be  rare  and  limited 
to  conditions  of  calm  sea  accompanied  by  light  winds.  Range  hole  phenomena 
were  not  detected  in  the  positioning  data  studied. 
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H.     CRITICAL  SIGNAL  THRESHOLD  (CST)  DERIVATION 

The  relatively  large  distances  and  varying  propagation  conditions  encount- 
ered in  the  Monterey  Bay  network  result  in  many  low  signal  strength  readings. 
It  is  important  that  a  method  be  found  for  using  these  readings  if  at  all  possible. 
Casey,  1982,  has  shown  that  low  signal  strength  ranges  can  be  used  if  corrected 
in  postprocessing  by  an  empirically  derived  relation  between  bias  error,  standard 
deviation,  and  signal  strength. 

The  recommended  procedure  to  establish  this  relation  is  to  perform  repeated 
range  measurements  over  a  known  length  baseline,  with  a  variable  attenuator  to 
artificially  vary  the  recei\ed  signal  strength.  These  attenuators  were  not 
available,  but  the  effect  was  achie\ed  by  ranging  over  relatively  long  baselines  of 
se\eral  different  lengths,  in  order  to  obtain  readings  over  a  spread  of  low  and 
high  signal  strengths.  Two  sets  of  range  readings  were  collected  for  different 
combinations  of  Master  and  reference  stations,  with  signal  strengths  varying  from 
about  5  to  25. 

For  each  Master  Reference  station  pair,  ranges  with  signal  strength  over  15 
were  a\eraged  to  estimate  a  mean  range,  free  of  low  SS  effects.  Differences  be- 
tween this  mean  range  and  all  low  SS  (  <  15)  ranges  were  computed,  representing 
the  additional  error  due  to  low  SS.  By  this  technique,  error  estimates  from  dif- 
ferent Master  Reference  station  pairs  could  be  combined  without  concern  about 
systematic  differences  between  equipment  pairs.  The  resulting  range  error  esti- 
mates were  grouped  by  signal  strength:  5-6,  7-8,  9-10  and  so  forth.  Mean  errors 
and  standard  dexiations  were  then  computed  for  each  signal  strength  group  as 
shown  in  Table  9.  As  expected,  both  the  mean  error  and  standard  dexiation  tend 
to  increase  as  the  signal  strength  decreases. 

The  two  data  sets  were  compared  to  ensure  repeatability  of  the  results  before 
deri\  ing  the  empirical  relation.  In  Table  9  we  can  see  the  results  of  the  statistical 
tests  for  the  two  data  sets.  They  are  treated  as  samples  from  one  population,  but 
each  with  their  own  estimated  means  and  standard  de\iations  and  examined  with 
the  Fisher  -  Behrens  statistical  test  [Hamilton,  1964]  at  a  95" o  confidence  lexel. 
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as  follows.  The  T  test  limit  was  taken  from  the  Student's  t  distribution  tabic  wilh 
degrees  of  freedom  calculated  by: 


dof  = 


2  2 


+ 


n\[i^  -  \)  nl[}i2-  \) 


The  calculated  value  for  the  test  is  gi\  en  by: 

(.?!  --T2)-  (^/,  -  A/2) 


u  = 


Where  a\  and  cr^  are  the  estimated  variances,  /j,  and  /^2  ^^e  the  estimated  mean 
errors,  and  n^  and  n^  are  the  number  samples  in  the  first  and  the  second  data  sets 
respectively. 

We  examine  the  hypothesis: 

Uo'-^x  —  Xj  ^  ^    The  t^^'<-'  ^''^'^'^  sets  belong  to  the  same  population. 

Against  the  allernatiNe 

//,:Jv,  —  ^2  ^  ^  The  two  data  sets  don't  belong  to  the  same  population.  Table 
9  shows  the  deri\  ed  statistics  from  the  two  data  sets  and  the  limit  of  the  T  test  in 
each  case.  E\ery  calculated  value,  for  the  Fisher-Behrens  statistical  icsl.  that  is 
bigger  than  the  T  \  alue  means  that  the  H^  hypothesis  is  rejected  so  something  else 
going  on  in  the  data  sets. 

In  the  column  'Calculated  value',  the  number  is  deri\ed  using  the  incans  and 
the  standard  deviation  of  the  measurements.  We  can  see  that  when  SS  is  '^-<^  and 
17-18,  the  calculated  value  is  bigger  than  the  T  test  limit,  which  means  that  in 
these  cases  the  null  hypothesis  is  rejected. 

The  number  gi\en  in  the  column  'Calculated  Milue  with  increased  cr  '  is  a 
combination  of  the  standard  deviation  of  the  long  term  variation  in  the  instru- 
ment and  its  observed  standard  deviation.    As  we  can  see  the  result  is  a  total 
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Table  9.       STATISTICAL  TESTS  FOR  THE  TWO  SAMPLE  DATA  SETS 


s.s 

1st  data  set 
(mean  error  /  a) 

2nd  data  set 
(mean  error/  o  ) 

T  test 
limit 

Calcu- 
lated 
value  (u) 

Calcu- 
lated 
value  with 
increased 

fX 

5-6 

15.1  4.5 

13.5-5.3 

1.98 

2.48 

1.78 

7-8 

8.1/3.7 

7.8  3.8 

1.98 

0.63 

0.41 

9-10 

4.5/3.0 

3.8/2.6 

1.98 

1.68 

0.98 

11-12 

1.9  1.7 

1.9/1.7 

1.98 

0.28 

0.19 

13-14 

0.3  1.8 

0.3,0.7 

1.98 

0.28 

0.2 

15-16 

-0.0  1.5 

0.1  0.6 

1.98 

0.98 

0.37 

17-18 

0.8  1.4 

0.3/0.7 

1.98 

3.35 

1.17 

statistical  agreement  of  the  data  sets.  We  conclude,  therefore,  that  they  belong  to 
the  same  population. 

The  next  step  is  the  derivation  of  the  equation  that  will  permit  the  use  of  the 
lower  SS  more  effecti\ely.   The  cc|uation  deri\ed  was 

y^(  71.093  _27.206) 

llLX 

Where  y  represents  the  mean  error  in  meters  and  x  the  low  signal  strength 
and  the  given  numbers  represent  meters.  The  resulted  y  \alue  represents  the  ob- 
served minus  the  calculated  (C-O)  range  correction. 

Figure  13  shows  this  resulting  best  fitted  cur\e.  The  correlation  for  the  case 
is  more  than  0.999  which  indicates  that  the  resulting  equation  is  \ery  close  to 
perfectly  modelling  the  data.  This  equation  can  be  included  in  the  software  as  a 
correction  for  the  lower  SS  rather  than  rejecting  them  below  the  CST. 
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Figure   11.      FllTED  Cl!R\  E  TO  THE  COINIBINED  DATA  FOR  lESTS  1  AND 


Additionally,  due  to  the  increased  standard  de\iation  of  the  ranges  at  lower 
SS  a  higher  value  for  standard  deviation  of  a  single  measurement  should  be  used 
in  the  weight  matr  x  formulation. 

I.     VARIATION  OF  THE  STANDARD  DEVIATION  WITH  TIME 

Eigth  sets  of  range  measurements  were  collected  over  the  known  distance  be- 
tween station  HA^l  S  and  geodetic  station  RANGE  7.  The  choosen  distance  was 
short  (2341   m)  to  eliminate  the  effect  on  distance  of  the  \ariation  of  weathci' 
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conditions.  These  8  sets  were  collected  on  8  different  days  and  they  have  length 
of  about  1  hour,  with  3  second  intervals  between  range  measurements.  Table  10 
gives  the  bias  of  these  measurements  with  respect  to  the  known  distance  and  the 
standard  deviation  of  a  single  measurement.  The  results  from  Table  10  show  that 
the  system  has  a  standard  deviation  which  is  indicative  of  its  precision  rather 
than  its  accuracy.  In  addition,  it  has  a  variable  bias.  The  total  standard  deviation 
for  the  system  must  reflect  both  these  components. 

Table  10.       STANDARD     DEVIATION     AND 
BIAS  OF  RANGE  ERROR 


NUM 

Standard  de\i- 
ation 

Bias 

1 

0.285 

-0.725 

0 

0.206 

1.447 

3 

0.255 

1.324 

4 

0.266 

1.531 

5 

0.374 

0.577 

6 

0.856 

0.007 

7 

0.576 

0.501 

S 

0.541 

1.173 

The  biases  shown  in  Table  10  were  used  statistically  to  show  that  the  stand- 
ard de\iation  of  one  measurement  is  2.54  meters,  which  concurs  with  the  value 
deri\ed  by  evaluation  of  g]  estimates  in  section  C  of  this  chapter.  The  procedure 
followed  for  this  derivation  was  a  simple  multiplication  of  the  standard  deviation 
of  the  mean  bias  found  from  the  data  set  by  the  square  root  of  the  number  of 
data  sets.  The  user  shouldbe  aware  that  the  standard  deviation  can  be  larger  on 
board  a  vessel  in  a  rough  sea  due  to  mast  movements  and  the  selcctivitv  of  the 
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stations  that  the  system  makes(non  simultaneous  measurements  to  the  reference 
stations).  The  bias  variation  from  day  to  day  makes  the  calibration  of  the 
stations  time  consuming  because  a  one  day  calibration  will  not  remove  all  the 
bias. 

Table  11.       MEAN     DIFFERENCES     AND     STANDARD     DEVIATIONS     FOR 
STATIONS  HAYS  AND  TREVOR 


DIS- 
TANCE 

Mean  differ- 
ence (Hays) 

Standard  devi- 
ation (Hays) 

Mean  differ- 
ence (Tre\or) 

Standard  de- 
viation 
(Trevor) 

SHORT 

0.55 

0.8 

- 

- 

ME- 
DIUM 

1.02 

0.55 

0.91 

0.48 

LONG 

1.27 

0.75 

1.48 

0.68 

V.LONG 

1.2S 

1.2 

0.84 

0.95 

J.     VARIATION  OF  THE  STANDARD  DEVIATION  WITH  DISTANCE 

The  second  thing  was  the  variation  of  the  standard  de\"iation  with  the  dis- 
tance or  with  S.S.  when  it  is  abo\e  the  CST.  Se\eral  problems  limited  the  data 
collection,  like  the  una\'ailability  of  an  attenuator,  the  poor  network  of  known 
positions  around  the  Bay  and  the  low  ele\ation  of  possible  calibration  sites.  The 
data  were  sufficient  only  for  the  e\'aluation  of  station  HA^'S  and  partially  of 
station  TREVOR. 

Table  1 1  shows  the  mean  difference  of  the  measured  distances  from  the  true 
distances  and  the  standard  de\iations  of  each  set  of  range  measurements.  We  can 
see  that  there  is  variation  of  the  mean  error  with  distance  in  both  stations  but  no 
indication  exists  that  a  function  can  be  deri\ed  because  the  mean  differences  with 
the  izi\en  standard  de\iations  are  not  statisticallv  different.  Actually,  a  functional 
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variation  of  the  standard  deviation  with  distance  was  not  expected,  when  oper- 
ating above  CST. 
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V.     DISCUSSION  AND  CONCLUSIONS 

The  Mini-ranger  networlc  in  the  Monterey  Bay  was  installed  for  research 
purposes  by  MBARI  and  NPS  and  was  examined  from  various  aspects.  Useful 
results  were  drawn  that  can  help  in  the  understanding  of  the  operation  of  trans- 
ponder type  systems  in  other  areas,  and  that  will  make  the  present  network  more 
useful  when  operating  at  low  SS. 

Figure  12  shows  the  o\erall  coverage  of  the  Mini-ranger  with  all  the  hmita- 
tions  taken  into  acount,  except  that  weather  conditions  can  attenuate  the  signal 
and  in  turn  reduce  the  maximum  range  of  the  stations.  For  the  maximum  range 
derivation  a  height  of  the  master  station  on  the  ship  was  assumed  to  be  10  meters. 

From  this  network  and  the  data  processed,  the  following  conclusions  can  be 
drawn: 

•  For  the  calibration  of  such  systems,  if  one  wants  the  maximum  accuracy  that 
can  be  pro\ided.  a  bias  taken  from  single  obser\'ations  over  known  distances 
or  with  the  baseline  method  may  be  not  enough,  it  is  better  to  use  repeated 
obser\"ations  and  a  mean  corrector  derivation. 

•  The  positions  of  the  stations  in  the  Bay  were,  with  the  exception  of  station 
HANK,  determined  to  a  sufficint  level  of  accuracy  such  that  they  were  not 
responsible  for  some  of  the  large  range  standard  deviations. 

•  The  software  of  the  user  is  important  in  the  position  derivation  and  can 
cause  unexpected  errors  if  not  properly  validated. 

•  The  meteorological  conditions  in  the  Bay  can  cause  attenuation  of  the  signal 
and.  in  turn,  reduce  the  maximum  range  ability.  The  range  error  from  the 
delay  of  the  signal  due  to  meteorological  conditions  does  not  exceed  the  1.2 
meters  in  the  extreme  conditions  and  is  smaller  than  0.2  meters  under  the 
usual  operational  conditions. 

•  Range  holes  in  the  Bay  can  occur,  but  should  be  rare  due  to  usual  windy 
Bay  conditions  and  the  modulated  surface  from  the  swell. 

•  We  have  shown  that  ranges  measured  at  low  SS  need  not  be  rejected,  but 
can  be  corrected  using  an  appropriate  error  function.  This  function  may  vary 
from  svstem  to  svstem. 

•  We  show  that  the  standard  deviation  can  be  considered  to  consist  oi^  two 
components,  the  unpredicted  bias,  which  is  smaller  when  the  system  is  well 
calibrated  and  the  measured  standard  deviation  of  the  data.  Generallv  the 
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Figure  12.      MINI-RANGER  COVERAGE  IN  MONTEREY  BAY 

a  of  one  measurement  under  ideal  conditions  is  a  little  greater  than  2  metcis. 
It  is  more  proper  to  use  (7  =  3  meters  when  the  system  operates  at  sea. 
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•    There  are  no  indications  for  variation  of  the  standard  deviation  with  dis- 
tance (or  SS)  above  CTL. 
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APPENDIX  A.     COMPUTER  PROGRAM  LISTING 

REM  PROGRAM  "LESS"  FILE:   NICKS. BAS 

'     NICK  KRIONERITIS      Date:  8-8-89 
Modified   10-17-89  by  KJS 

1  ^^  .' .  y ^  y « ..'^  .u  ^.  ^.  y .  y^  .t^  y^  y^  y ^  y ^  y .  y «  ^*.  .'^  y ^  y ^  ^  .t^  ^,  ^.  y^  ^  .«,  y^  y^  ^.  y .  ^.  y^  y^  ^u  y^  ^^  y^  ^^  y .  y^  y^  yu  y^  y^  y .  y .  y .  y ^  y .  y. 

REM  This  program  calculates  the  optimum  position  of 
a  ship  using  distances  from  known  stations.  It 
serves  the  solution  for  a  miniranger  network  with 
accuracy  better  than  10  centimeter. 

t  V^  Vr  Vr  Vr  Vr  ?V  V'  Vc  Vc  Vr  ^'  Vr  Vr  V'  Vr  ~V  Vc  Vc  Vr  Vr  Vr  Vr  Vr  Vr  *V  V'  Vr  ^  ~V  V'  ^  V'  Vr  ^'  "^r  V'  Vr  V'  Vr  Vr  V'  V*  Vc  Vr  Vc  V'  'V  Vc  V'  ^  'V  ^V  V'  Vr  Vr  'V  Vr  Vr 

REM  The  subrutine  SOLVERANGX  was  changed  from  CDR  K. SCHNEBELE 
program  "HYDROPLOT". 

f  -J*  y .  ^.  y .  y.  .J^  ^'.  .1^  ^f^  ..f .  ..t.  J^  .J..  .1.  .'.  J.  ^1.  J>~  .*.  .u  .*..  .J..  k<-  .J.  J«  .J.  .'.  .^,  ^.  »t.  .f.  ^t.  ^.  ^.  j^  j^  ^.  ^.  ^^  J.  ^.  ^.  ^«  y«  ^.  ^«  ^^  ^.  «*.  ^1,  y.  ^«  y.  ^.  ,j.  ,j.  .1^  ^<. 

$STACK  9216 

$ INCLUDE  "SOLRANGX.  SUB" 

$INCLUDE  "NICKMR. SUB" 

ON  ERROR  GOTO  ERRORHANDLER 

Print  "Enter  the  number  of  stations" 

Input  "N°o=",N°o 

Lprint  "Number  of  stations="  ,N°o 

DECLARATION  OF  THE  VARIABLES 

AP/ZCl):  Approximate  position  of  theship  in  meters  (UTM) 

SP//(2):  Station  positions  in  2D  array  in  meters  (UTM) 

Rj-'/Cl):  Ranges  observed  on  the  receiver  in  meters 

LX#(1):  Stations   Eastings  for  the  scale  factor  affection 

X.:/(l):  False  Eastings  of  the  stations 

Y?K1):  Northings   of  the  stations 

Ziril):    Elevation  s  of  the  stations  from  MSL 

SF#(1):  Scale  factor  calculated  for  each  case 

W!(l):  Diagonal  Weight  matrix  assuming  indepedence 

among  the  measured  ranges 
C0DE%(1):  Index  numbers  of  the  4  stations  in  each  record 
AllSta/K2):  UTM  Coords  (x,y,z)  of  6  stations  possible  to  use 

Vectors  for  S0LVERANG2  but  not  used  in  this  version 
Cm0s!(l):  Bias  corrections  for  each  of  6  stations 
Wgts!(l):  Weight   factor  for  each  of  6  stations 


DIM  AP/Zf  1:  3) 

DIM  DYNAMIC    SP#(  1:  N%,  1:  3)    ,    R//(  1:  N?o)    ,    LX#(1:N%) 

DIM  DYNAMIC   X#(l:N?o)       ,    Y//(  1:  N%)    ,    Z//(1:N%) 

DIM  DYNAMIC   M,-'/(l:N°;)       ,    SF//(  1:  N°i)    ,    BRes#(l:N%)    ,    W!  ( 1:  N%) 

DIM   code%(l:4),AllSta#(l:  6,1:  3),Cm0s!  (1:  6),Wgts!  (1:  6) 
t 

PRINT  "FREE  ARRAY  &  STACK  AT  STARTUP",  FRE(-l),  FRE(-2) 

CALL  GetStations(AllSta#(),CmOs! (),Wgts! ()) 

Lprint  CmOs! ( 1); CmOs! ( 2); CmOs! ( 3); CmOs! (4); "  CmOs" 

CALL  SetFiles 

Print  "Enter  assumed  position" 

Input  "AX=",AP#(1) 
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Input  "AY=",AP//(2) 
Input  "AZ=",AP//(3) 

LPRINT  "Assumed  Position"; 

For  i%=l  to  3 

Print  "Assumed  position  in  UTM  (X,Y,Z)="  AP#(i%) 

LPRINT  USING  '' iHHHHHHHh  U'' ;    AP#(i%); 

Next  i% 

REM  Aproximations  for  calling  the  subrutine  SOLVERANGZ'"^ 

SFAC!=1.0       'Approximations  for  the  subrutine  for  the 

SHGTy/=AP//(3)     '     corrections  that  have  already  done 

AP#(3)=0.dO 

LM7'/=1.  1         'Convergence  Limit  test  (meters) 

Rm,'/=6368459  'Mean  radius  of  Earth  in  meters     '"^  0.9996 

REM: Declaration  of  Weight  matrix  * 

For  i%=l  to  N% 
W!  (i%)=.  25 


'Diagonal  matrix  assuming 
' independent  measurements 

Next  i% 

Beep  2:  CLS 

PRINT  "Put  Printer  online  in  COMPRESSED  PRINT  Mode!" 
PRINT  "   Hit  any  key  when  ready  to  continue" 
WHILE  NOT  INSTAT 
LOCATE  2,37 
PRINT  "_" 
WEND 
Print  INKEY$ 

'   LOOP  TO  READ  AND  PROCESS  ALL  DATA  UNTIL  END  OF  FILE  HIT 

'WHILE  NOT  EOF(l) 
For  k%=l  to  500 

CALL  READNICKMR(sec<Sc,code%( )  Mi )  ,X//( )  ,¥//( )  ,ZiK  )  ,AllSta#( ) ) 

'  it'k-k'k'!-'k-k'k-!!:-k'k-k'k-';-it-!i:-!<:-k-k'k-!!:i--k-k-k--k'k'k'k'k-!:'!-'k'l"k 

REM  Correction  for  slope  distance  - 

For  i%=l  to  N% 

If   {MUX)    >    l.dO),    Then 
DHGT//=Z//(i%)-SHGT// 

HGT#=R,^i(i%)  2-DHGT//   2 

IF   HGT7,'<   O.dO,    Then 

PRINT  "Warning   -   at   time   "; sec&; "   DHGT  exceeds   Range" 
PRINT  "   Hit   any  key  to   continue" 
While   Not    INSTAT 
Delay   0.  2 

n       .      ^      IT      II 

Print      .     ; 
Wend 
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Print  INKEY$ 
Else 

R^/(i%)=SQR(HGT#) 

End  If 
Else 

R#(i°o)  =0.dO 
End  If 
Print  "Ranges  corrected  for  slope="  R//(  i%) 
Next  i% 

REM  Calculation  of  the  scale  factors  "''=■ 

AXP#=AP#(1) -500000 
For  i%=l  to  N°o 

LX#( i%)=X#( i%) -500000 

M#(i%)=AXP#  2  +  AXP//*LX^/(i%)  +  LX/ZCi^)  2 

STiKi%)=0.9996  +   .999 6'>M# (  i% )  / (  6'- Rm-'Z ) 

Print  "Scale  Factors  for  UTM="  SF^/(i?o) 

'   Lprint  "Scale  factors  for  UTM="SF//(  i%) 

REM  Correction  for  the  scale  factors  -'" 

R#(i%)=R#(i%)^'^SF?;i(i%) 

Print  "Corrected  ranges="  R#(i%) 

Lprint  "Corrected  ranges="  R^r'K  i/o) 
Next  i% 

REM  Formation  of  station  position  matrix   "''' 

I      .J.  ^.  .f.  -'r  Vr  V'  *"''  "*''  ^'  "''  ^*'  '*'  "''  -^"^"^'^^  Vr  -'-  o^'  -'^'  ^'  -V*^'  ^*'  *'''  ^*'  ^-  ^'  "*'-*'  ^-  ^*'  -V  ^*'  ^*'  ^*'  '^'  ^''  ^V  "^r  "V  "  V 

For   i%=l   to  N?b 

SP#(i°o,l)=X#(i%) 

P,-'/(i?i,2)=Y,^r'(i%) 

SP//(i%,3)=0.dO 

Next    i% 
t 

CALL  SolveRangX(AP#( )  jSP^/C )  ,R,^/( )  ,CmOs!  ( )  ,N%,SFAC!  ,   W!  ( ) ,    LM//    ,BRes-:'/( ) ,. 

sxx!  ,    syy!  ,    sxy!  ,    variance!  ) 
CALL  RiteSoKsecajAP/iO  ,code°o()  jR/ZO  ,BRes#()  ,sxx!  ,syy!  ,sxy!  , variance!  ) 

LPRINT  USING   "   iHHHm  imiHHHh  il^  ##M#/A^/.#";    sec&   AP//(1)    AP#(2); 

For   i'o  =   1   to  4 

LPRINT  USING   "   ^H^  iHHHh  i^  HHhir;    code%(i?i)   R#(  i?o)    BRes//(i%); 

Next    i° 

LPRINT  USING  "  im.UU-i^   +###.#  UU-i^;    sxx!  syy!  sxy!  variance! 

'WEND 
t 

Next  k% 

CLOSE  y/1 

CLOSE  #2 

END 

ERRORHANDLER: 

PRINT  "FREE  AT  ERROR" ,FRE( -1) ,  FRE(-2) 

PRINT  "AN  ERROR  TYPE",  ERR,  "HAS  HAPPEN  AT",  ERADR 

END 
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SINCLUDE  "openfile.  sub" 

*-jV-!V-v-.v-;oViV-,v-iV,v-;"V->v  GetStations  Subroutine  *********iWov  kjs:Oct  89  -"=■ 

file: NICKMR.  SUB 
SUB  GetStations(AllSta#(2),CmOs!  (l),Wgts!  (1)) 
LOCAL  i?4,  s$,  code% 

CLS:  PRINT  "insert  Disk  with  File  ST. DAT  in  Drive  B: " 
Print  "  Hit  any  key  to  continue  when  ready"; 
While  Not  INSTAT 
Locate  2,37 
rmt   ! 
Wend 

Print  INKEY$ 

CALL  0PENFILE("I",1,"B:  ST.  DAT", 128) 
For  i%   =  1  to  6 
INPUT  in,    code%,AllSta#(i%,l),AllSta#(i%,2),AllSta7'/(i%,3),CniOs!  (i%),_ 
Wgts!  (i?o),s$ 
Next  i% 
CLOSE  in 
END  SUB 

Vr,vvrVc-.v--v--v,v,v--vvr  SetFiles  Subroutine  *-v,v,v,v.v,v.v,v.v.v-'«v.v-vVfVf,vvr  ]^js:    Oct  89  'V-.v-v-v-.v-v-;.- 

f  ile:  NICKMR.  SUB 

opens  input  data  file  of  Miniranger  values(as  #1),  and  opens  output  file 
for  the  computed  results  (as  ,-'/2) 

SUB  SETFILES 

LOCAL  s$,txt$,FiName$ 

CLS:  s$  =  "  (Drive:  path  name,  ext) 

PRINT  "Enter  MiniRanger  data  file  for  input"  +  s$; 

INPUT  FiName$ 

CALL  0PENFILE("l",l,FiName$,128) 
t 

PRINT  "Enter  solutions  file  for  output"  +  s$; 
INPUT  txt$ 

CALL  0PENFILE("0",2,txt$,128) 
PRINT  ?'/2,  "Solutions  from  file  ";  FiName$ 
END  SUB 
Vrvrvrvrvrvrvrvrv-V'vr  ReadNickMR  Subroutine  VrVr^Vr.vvrVcv-VfVr:?cv-Vf-v,vvrVr.v  kjs:  Oct  89  -'-.v-vvr.vvr 

file:  NICKMR. SUB 
specific  to  read  mini-ranger  files  (from  #1)  from  Krionertis  thesis  work 
with  time  in  sees  &  4  station  numbers  &.   4  ranges. 

SUB  READNICKMR(sec&,code°o(  1)  ,rang5'/(  1)  ,X//(1)  ,Y//(  1)  ,Z#(1)  ,  AllStayi(  2)  ) 
LOCAL  s$,i%,nucode°o( ) 
DIM  DYNAMIC  nucode?o(l:  4) 

LINE  INPUT  in,    s$ 
sec&=VAL(MID$(s$,8,7)) 
For  i%  =  1  to  4 
nucode°o(i°o)=VAL(MID$(s$,ll+4''-i%,4)) 

rang#(i%)=VAL(MID$(s$,22+9'H%,9)) 
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'check  nucode?i( ) ,  if  not  same  as  code%( )  then  get  new  station  coordinates 
I  for  the  X//(),Y//(),Zy/()  arrays 

If  (code%(i%)  <>  nucode%(i%)),  Then 
code%(i%)  =  nucode%(i?0 
X//(i%)  =  AllSta//(code%(i%),l) 
Yy/(i%)  =  AllSta,-'/(code%(i%),2) 
ZiKiX)   =   AllSta#(code%(i%),3) 
End  If 
Next  i?^ 
END  SUB 
'f-;-iVv--v-ViVVr^v**Rit:eSol  Subrutine  ****-'-**->'^**********---**-'^**->''->'^   kjs:  Oct  89-'^"""'^ 

File:  NICKMR.  SUB 

Writes  Miniranger  Solution  and  error  statistics  to  File  #2 

SUB  RiteSol(sec6c,AP#(l),code%(l),R//(l),BRes#(l),sxx!  ,syy!  ,sxy!  ,vari!  ) 

LOCAL  i% 
t 

PRINT  n,   USING  ''Umm  iHHHHHHhi^  iHHHmihir;    secSc  AP//(1)  AP#(2); 
f$  =  "#//  UUU.i^   +M.#" 
For  i?;  =  1  to  4 

PRINT  y/2,  USING  f$;  code%(  i%)  R#(i%)  BRes//(i%); 
Next  i.% 
PRINT  n,    USING  "  m.U  ihhU   +/;//.##  UU.W;    sxx!  syy!  sxy!  vari! 
END  SUB 


File  =  SOLRANGX.  SUB 

vr  v.-  Vf  -V  vr  i:   SOLVE  RANGE  POSITION  (ver  X)  '-   ^< 

Returns  grid  x,y  in  first  2  elements  of  apsn7"'/( ) 


SUBROUTINE 


INPUTS: 

apsn7^K3) 

spsns?r 

(nsta°o, 
ORSrKnst 
rangcor! 
nsta°o 
scalefac 
wgts! (ns 
dx  1  im// 

OUTPUTS: 
apsn//(3) 

BRes//(ns 
sxx! 

syy! 
sxy! 

variance 


3) 

a?o)  = 
(nst 

!   = 

ta%) 


approximate  position  as  grid  x,y,z 
station  positions  in  2D  array  as  grid  x,y,z 
(array  has  nsta°o  rows  by  3  columns  for  x,y,z) 
observed  ranges  to  each  station  in  spsns#  row  order 

a°o)  =  C-O  corrections  for  each  range  reading 
number  of  stations/ranges 
grid  scale  factor  (grid: true  ratio) 
=  weights  for  each  rang^  (1/sigma 

convergence  limit  to  stop  iterating  position 
solution  (given  in  grid  units  -  e.g.  meters) 


=  x,y  position  returned  in  1st  two  elements 
(elevation  of  antenna,  z,  unchanged) 
ta?o)  =  residuals  on  observed  ranges  (ob-computed) 
=  variance  in  x  (from  Inv(ATWA)) 
=  variance  in  y 
=  covariance  of  x  &  y 
!  =  est  of  sigma-o  (2?i\^]BXiWi^tE3-2) ) ) 


$INCLUDE  "LUDCMP. BAS" 
$INCLUDE  "LUBKSB.BAS" 
SUB  So  lveRangX(apsn#(l),spsns/K  2  ),ORS,-'/(l),  rangcor!  ( 1)  ,nsta%, scalefac!  , 
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wgts!  C 1)  ,dxlim//,BRes#(  1)  ,sxx!  ,syy!  ,sxy!  , variance!  ) 
LOCAL  A#( )  ,spsn#( )  ,ORa//, Ax//, Ay#,Bi#,ATW#( )  ,ATWB#( )  .ATWA/ZC )  ,Indx%( ) , 

suiTi#,iteration%,i%,  j%,k%,YCol#( )  ,ATWAInv#( )  ,xdx%( ) ,    Btemp/ZC) 
DIM  spsn//(l:3),    ATWB#(1:2),    ATWA//(  1:  2, 1:  2) ,    Indx%(  1:  2)  ,xdx%(  1:  4) ,_ 
Btemp,-'/(l:4) 

i%=l 

For   j%  =   1   to  nsta?^ 

If   ORs//(j%)    >    l.dO,    then 

xdx?o(  i%)   =  j/o  'xref  to  non-zero   ranges    (range  >    Im) 

INCH   iX 
End   If 
Next   j% 
nranges/o  =   i%-l 

DIM  A,-'/(l:nranges%,l:2),    ATW//(  1:  2, 1:  nranges%) 

SolveAgain: 

For    i%  =    1    to  nranges'/o 

ORa#=ORS#(xdx%(i%))   +  rangcor! (xdx%( i%) ) 
For   2%  =   1   to   3 

spsn,-Kj°o)   =  spsns//(xdx%(i%),j%) 
Next   j/o 

CALL  RangeLOP(  apsn#( ) , spsniK  ) ,  scalef ac!  , ORa// ,  Ax// ,  Ay// , Bi//) 
A//(i%,l)=Ax// 
A//(i%,2)=Ay# 

BRes//(  i/o)=Bi//  'Residual   (o-c)    on  each   range 

Next    i°o 

'      Compute  ATW   02  ,nranges?;|  ,    &  ATWB   02  | 

For   j%  =   1   to   2 
suni//=0.  do 
For   i/o  =   1   to  nranges?b 

ATO/Kj%,i°o)   =  A//(i?o,j%)"Wgts!  (xdx%(i%)) 
sum//  =  ATW7/(  j?i ,  i%)-'fBRes//(  i°o)+suni// 
Next    i% 
AT\N'B//(j?o)=sura// 
Next   j°o 

'      Compute  ATWA  02,2 | 
t 

For   j°o  =    1   to   2 
For    i°o  =    1    to   2 
sum//  =   0.  do 
For   k%  =    1    to  nranges% 

sum^/  =  ATC//(j%,k°o)''-A//(k%,i%)+sum// 
Next   k% 
ATCA//(j°o,i%)=sum// 
Next   i?i,j°b 

CALL  LUDCMPC  ATU"A//( ) ,  2  ,  Indx%(  ) ,  i%) 
CALL  LUBKSB(ATOA//( )  ,2,Indx°o( )  ,ATOB//( )) 
sum//=0.  do 
For   i%  =1   to   2 
apsn//(i%)=apsn/i(i?;)+A'nvB//(i%)         'shifts   position  by  dx,dy  which  has 
sura//  =  ATWE/ECi-^o^um//  'been  returned   as   AT\\'B(,)    from  LUBKSB 
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Next  i% 
sum#  =  SQR(sum//) 
If  sum//  >  dxlim//,  Then 
INCR  iteration?^ 
If  iteration%  >  7,  Then 
Beep  3 

Print  "Position  did  not  converge  --NO  SOLUTION!  "; 
Print  Using  "Delta  Position  =  +#.  ####meters";  sum# 
sxx!  =  0.0:  syy!  =  0.0:  sxy!  =  0.0:  var!  =  LOG10(sura#) 
EXIT  SUB 
End  If 
GoTo  SolveAgain 
End  If 

Compute  variance-covariance  values  from  inverse  of  ATWA 
sxx!  =  sigma-x  2 

syy!  =  sigma-y  2 

sxy!  =  sigma-xy 

ERASE  ATW//,  spsn//  'make  some  room  in  memory 

DIM  YCol,-'/(l:2),  ATU'AInv#(l:  2, 1:  2) 
For  j%  =  1  to  2 

For  i%  =  1  to  2  'identity  vector,  YCol,  to  find 

YCoiy/Ci^O  =  0.  do  'inverse  using  back-substitution 

Next  i% 
YCol#(j°c)  =  l.dO 

CALL  LUBKSBC  ATWA/ZC ) ,  2  ,  Indx%( ) ,  YCoiy/( ) ) 
For  i?;  =1  to  2 

ATU^AInv//(i?^,j%)=YCol#(i%) 
Next  i% 
Next  j% 
sxx!  =  CSNGCATWAInv/ZCljl)) 
syy!  =  CSNG(AT\vAInv,-'/(2,2)) 
sxy!  =  CSNG(AT\v'AInv#(l,2)) 

Compute  estimate  of  unit  v^^igtix^/iEBaiS^'^ftBS 
Note  that  range  residuals  are  in  BResT'/() 

If  nranges°i,  =   2,    Then 

variance!  =  9999  'signifies  no  variance  computable 

Else 

sum//  =0.  do 

For   i%  =   1   to  nranges°o 

sum//  =  Bres//(i°o)  2"Wgts!  (xdx%(  i%) )   +  sum// 

Next    i% 

variance!  =  CSNG(sum///(nranges%-2)) 
End  If 

If  nranges?o  <  nsta°o,    Then  'Set   residuals    into  proper   channel 

For   i°o  =   1   to  nsta°b 
Btemp//(i%)   =  0.  dO 
Next    i% 
For   j%  =   1   to  nranges% 

Btemp//(xdx%(j%))   =  Bres//(j%) 
Next   2% 
For   i°o  =   1   to  nsta?o 
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■jV-j'oV 


If  nranges%  =  2,   Then  Bres#(i%)   =  0. dO 

BresiKiX)   =  Btemp#(i%) 

Next   i% 

End   If 
t 

END  SUB 
*  *   ,v  *  vr  Vf  RANGELOP  *  *  SUBROUTINE  *  *  *  *  Vc  -v  *   *  *  Vr 

Computes  A  matrix  components  Ax  &  Ay  and  grid  range 
for  input  to  multiple-LOP  solution  algorithm. 

Inputs:  apsn?'/( )  =  approx  position  x,y,z  coordinates  (grid) 
spsn,"'/( )  =  station  x,y,z  coordinates 
scalefac! =  grid  scale  factor  (point  or  line) 
ORay/  =  observed  slant  range 

Output:  Ax#  &.   Ay//  such  that  ( Ax)dx+( Ay)dy  =  OR-CR 

OraCRT'/  =  Observed  minus  Computed  range  at  apsnT*/ 
Note:  Observed  range  (OR)  has  been  corrected  for  elevation 
differences  and  grid  scale  before  comparing  with  CR. 

SUB  RangeLOP(apsn?'K  l),spsn/Kl)  , scalefac!  ,_ 

ORa//,Ax//,Ay//,OmCR?/) 
LOCAL  i%,CR#,dx#() 
DIM  dx#(l:  3) 

For  i%   =  1  to  2 

dx//(i%)  =  apsn//(i%)-spsn,-'/(i%) 
Next  i% 

CRy/  =  SQR(dx//(l)  2+dx#(2)    2) 

Ax#  =  dx#(l)/CR,-'/ 
Ay?/  =  dx//(2)/CR# 

dx//(3)=0.d0 
OmCR//  =  ORar-'^scalefac!     -   CR# 
END   SUB 


file:OPENFILE.  SUB 

Opens  a  disk  file  using  TBasic  commands.   Preforms  error  checking  and 
allows  recovery  from  common  mistakes. 

INPUTS:   modes$  =  One  character  designation  for  filemode 

"O"  =  for  sequential  output  (to  be  written  to) 
"l"  =  for  sequential  input  (to  be  read  from) 
Csee  Turbo  Basic  Ref,  pg  282,  for  other  choices | 

FileNumX  =  file  number  to  be  used  in  read/write  calls 
FileName$  =  Drive: /path/name. ext  of  file 
RecLength?i  =  record  size  in  bytes,  1  to  32767,  typical  =  128 

SUB  OpenFile(mode$,FileNum°o,FileName$,RecLength%) 

LOCAL  errnum°o,  ans$ 
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OpenSequence: 

ON  ERROR  GoTo  ErrorHandling 

OPEN  mode? ,FileNum%,FileName$ ,RecLength% 

ON  ERROR  GoTo  0 

EXIT  SUB 

ErrorHandling: 

Beep  2 

errnumX  =  ERR 
Select  Case  errnum°o 
Case  53,75,76 

Print  "FILE  or  DRIVE  PATH  Not  Found  for  . .  .  ",  FileName$ 

Input  "Re-Enter  Drive:  path  name. ext  >  ";  FileName$ 

Case  68,71 

Print  "DISK  DRIVE  NOT  READY  for  ",  FileNatne$ 
Case  Else 

Print  "Opening  File  Gave  Error  Number  ",  errnum% 
End  Select 
Input  "(R)etry  or  (H)alt  Program  ";  ans$ 
If  Ucase$(ans$)  =  "H" ,  Then 
Resume  HaltSequence 
Else 

Resume  OpenSequence 
End  If 


I 


HaltSequence:      'terminates  main  program 

Print  "Program  Halted!  !  !  " 

Print  "Subroutine  OpenFile  with  Error  Number  ",errnum?4 
Print  "   encountered  at  program  address  ",  ERADR 
END 
END  SUB 

' LU  Back  Substitution 

Ref:  Numerical  Recipies,  by  WH  Press,  page  37 

SUB   LUBKSB(  AZ/C  2)  ,N?o,  INDX%(  1)  ,B#(1)  ) 
LOCAL   ll%,i%,j?i,,sum# 

ii%=0 

For  i%=l  to  N% 
ll%=Indx%(i%) 

sum#=B#(ll%) 
B#(ll%)=B//(i%) 
If  ii%<>0,  then 

For  j/o=ii%  to  i%-l 

s  um#=s um#  -  A#  (  i% ,  j  % )  '■  B//  (  j  % ) 
Next   j% 
Elself   sum?'/<>0.  do   Then 

±  ±  10 —  X  /o 

End   If 

B#(i%)=sum-:'/ 
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Next  i% 

For  i%=N%  to  1  Step  -1 
sum#=B#(i%) 
If  i%<N%,  Then 

For  j%=i%+l  to  N% 

s um5'/=s um//  -  A#  (  i% ,  j  % )  *B# (  j  %  ) 
Next  j% 
End  If 

B#(i%)=sum?'//A#(i%,i%) 
Next  i% 
END  SUB 

' LU  Decomposition  Subroutine 

Ref:  Numerical  Receipes,  by  WH  Press,  pages  34-36 

SUB   LUDCMP(A#(2),N?i,Indx?i(l),D%) 

LOCAL  vv#( )  jDum//,  i%,  j%,k%,  imax%,aamax#,sum# 

STATIC  nmax?o,    tiny// 

nmax°i=10 
tiny#=l.  OD-20 

DIM  vv//(l:nmax%) 

D°o=l 

For   i%=l   to  N?; 

aamax//=0 

For   j°o=l   to  N% 

If  ABS(A/Ki%,j°o))>aamax//,   Then   aamax#=ABS(A//(  i/o,  j%) ) 

Next   j% 

If  aamax//=0,  Then 

Print  "Matrix  has  zero  row  it   "  i% 

Exit  Sub 

End  If 

vv#(i%)=  l/aamax#         'save  row  scaling  factor 

Next  i7o 
t 

For  j%=l  to  N% 
If  j/o>l.  Then 
For  i/o=l  to  j%-l 
Sum,^^=A?'/(i%,j%) 
If  i^>l  Then 
For  k%=l  to  i%-l 

Sumy/=Sum#-A//(  i%,k%)'"A#(k%,  j%) 
Next  k% 
A#(i%,j%)=Sum// 
End  If 
Next  i% 
End  If 

continuing  j  loop,  searching  for  largest  pivot 
aamax#=0 

For  i°c=  j%  to  N% 
Sum#=A#(i°o,j?o) 
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If   j°o>l,    Then 

For  k%=l   to   j%-l 

Suin//=Sum//  -  A#(  i% ,  k%)*A//(  k% ,  j%) 
Next  k% 
A#(i%,j%)=Suiny/ 
End   If 

Dum//=vv#(i%)*ABS(Sumj'^)         'figure  of  merit   for  pivot 
If  Dumy/>aamax//,    Then 
imax/o=i% 
aamaxy/=Dum# 
End   If 
Next   i% 


If   j%<>imax%,    Then 
For  k%=l   to  N% 

Dum#=A#(imax%,k%) 
A#(imax%,k%)=A//(j%,k%) 
A#(j%,k?0   =  Durn# 
Next   k?o 
D%=-D% 

vv?'/(imax?i)=vv//(  j%) 
End  If 


' interchange  rows  if  reqd 


'change  parity  of  D% 


Indx?^(  j/o)=iraax?o  'divide  by  pivot 

If  (j%<>N%),  Then 

If  A?'/(j%,j%)=  0,  Then  A//(  j%,  j%)=tiny# 

Diim#=l/A?'/(j%,j%) 

For  i%=j?^+l  to  N% 

A#  (  i% ,  j  % )  =A#  (  i% ,  j  %  )  '>Dum# 
Next  i% 
End  If 
Next  ]% 

If  A//(N?o,N%)=0,    Then  A,-'/(N%,N?0=tiny// 
END   SUB 


vfVfirVr^rVc-.r'j'rr- 


r  Vr  -,'{■  ■)'?  V-  Vf  Vr  Vr  -V  -V  -k  -  V 


1,   585260.161,  4092490.284,   51.7,  0.0,  0.25,  TREVOR 


2, 
3, 
4, 
5, 
6, 


6029^ 


4085231.399,   23.1,  0.0,  0.25,  WATS 


609863. 128,  4076611. 345, 

607621. 289,  4055915. 264, 

600434. 009,  4051260. 014, 

596669.490,  4051826.411, 


33.  9,  0.  0,  0.  25,  PACKARD 

137.  2,  0.  0,  0.  25,  HAYES 

10.  0,  0.  0,  0.  25,  DOPPLER 

134.  5,  0.  0,  0.  25,  HANK 
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APPENDIX  B.        GEOMETRICAL  ACCURACIES  OF  THE  MONTEREY 

BAY  MINI-RANGER  NETWORK 


HAYS 


Monterey  Bay  Vicinity 


Figure   13.      SIAIION  tOMDINAIlON  1  REVOR-WA  IS-IIAYS-HANK 
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HAYS 


Monterey  Bay  Vicinity 


Figure  14.  STATION 

TREVOR-PACKARD-DOPrLER-HANK 


COMDINATION 
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Figure  15.      STATION  COiMBINATION  WAIS-PACKARD-HAVS-HANK 
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Figure   16.      STAIION  COMBINATION  PACKARD-HAVS-DOPPLER-HANK 


61 


HAYS 


MoDlercy  Pay  Vicintty 


Figure  17.      STATION  COMBINATION  TRE\  OR-WATS-HANK 
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Figure   18.      SIAIION  COMBINAIIONTREVOR-PACKARD-IIANK 
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Figure  19.      STATION  COMBINATION  WATS-HAVS-HANK 
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Figure  20.      STATION  COMBINATION  TREVOR-PACKARD-HAVS 
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